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INTRODUCTION 


7L> 

This  is  the  final  report  for  ONR-Contract  N00014-J5-C-^®*i',  "Investi- 
gations on  the  Dynamics  of  MHD  Flows."  The  contract  has  been  concerned 
with  MHD  flows  in  diffusers  and  liquid-metal  brushes,  and,  to  a minor 
extent,  with  the  spectroscopy  of  Helium  plasma  flows  in  homogeneous 
magnetic  fields. 

In  Chapter  II,  an  exact  similarity  solution  is  presented  for  MHD 

flows  in  plane  diffusers  when  the  external  magnetic  field  is  homogeneous 

and  parallel  to  the  apse  line  of  the  diffuser  walls.  For  this  geometry, 

the  Induced  magnetic  field  is  in  the  axial  direction,  too,  so  that  the 

Lorentz  force  is  irrotatlonal,  V x x B)  = ^.  It  is  shown  that  the 

magnetic  field  changes  the  pressure  distribution  in  the  flow,  but  does 

not  directly  affect  the  velocity  field  of  the  flow  and  the  onset  of  flow 

separation.  MHD  diffusers  with  axial,  homogeneous  magnetic  fields  are, 

therefore,  not  of  advantage  for  MHD  generators  since  they  do  not  permit 

to  control  flow  separation  by  means  of  the  applied  magnetic  field. 

In  Chapter  III,  a similarity  theory  for  MHD  diffuser  flows  with 

viscous  stress  relaxation  is  presented  when  the  external  magnetic  field 

is  azimuthal,  B (r)  « r e-  (induced  magnetic  fields  are  neglected 
o w 

assuming  small  magnetic  Reynolds  numbers,  1) » together  with  a 

kinetic  discussion  of  viscous  stress  relaxation.  We  show  that  viscous 

stress  relaxation  changes  the  form  of  the  velocity  field  at  low  Reynolds 

numbers  R and  reduces  the  critical  duct  angle  6®  for  flow  separation  for 
2 

Hartmann  numbers  H £ 2R/3,  and  that  flow  separation  is  inhibited  for 
2 

H > 2R/3  at  any  duct  angle  0^  < it. 

In  Chapter  IV,  an  exact  similarity  solution  for  MHD  diffuser  flows 
(without  viscous  stress  relaxation)  in  an  azimuthal  magnetic  field 

I 


it 


r' 


i 


B (r)  “ r is  presented,  under  consideration  of  the  Induced  mag.ietlc 
o 0 

field  which  is  radial.  Since  the  magnetic  Reynolds  number  of  MHD  flows 
is  generally  not  large,  ' 1»  we  could  demonstrate  quantitatively  that 
the  effect  of  the  induced  magnetic  fields  on  the  velocity  distribution 
and  the  onset  of  separation  is  small.  This  result  justifies  the  neglec- 
tion  of  induced  magnetic  fields  in  various  previous  investigations  on 
MHD  flows. 

In  Chapter  V,  the  Stark-effect  in  Helium  plasmas  with  homogeneous 
magnetic  fields  is  evaluated  for  absorbtlon  spectroscopy  in  presence  of 
laser  radiation  of  advanced  generator  flows.  It  is  shown  that  the  laser 
radiation  changes  qualitatively  and  quantitatively  the  satellite  spectra, 
and  that  the  investigations  on  this  subject  in  the  literature  are  insuf- 
ficient. 

Considerable  analytical  and  numerical  effort  was  expanded  in  studying 
the  influence  of  the  Hall -effect  on  MHD  diffuser  flows.  Unfortunately, 
this  investigation  could  not  be  completed  in  time  for  this  final  report. 

An  account  of  the  latter  work  will  be  communicated  as  soon  as  the  doctoral 
research  of  S.  H.  Choi  has  been  completed. 


■ 
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NONLINEAR  BOUNDARY-VALUE  PROBLEM  FOR  SELFSIMILAR  MAGNETOHYDRODYNAMIC  DIFFUSER  FLOWS 


, ACROSS  HOMOGENEOUS  MAGNETIC  FIELDS* 

By 

H.  E.  Wilhelm  and  S.  H.  Choi 
Abstract 

Similarity  transformations  and  selfsimilar  solutions  are  derived 
for  the  incompressible  radial  flow  of  conducting,  viscous  fluids  across 
an  external,  homogeneous  magnetic  field  (parallel  to  the  apse-line)  in 
a diffuser  with  electrodes  in  the  planes  0 = ±6^  (cylindrical  coordi- 
nate system).  The  conducting  flow  across  the  (axial)  external  and 
Induced  magnetic  fields  Induces  radial  and  azimuthal  current  densities, 
of  which  only  the  azimuthal  current  density  produces  a net  current  flux 
I 0 through  the  electrodes,  presumed  that  these  are  connected  by  an 
external  circuit.  The  eigenvalue  problems  for  the  radial  velocity  and 
Induced  magnetic  field  amplitudes  f(e)  and  ij;(0)  are  solved  in  closed 
form.  The  external  and  induced  magnetic  fields  are  shown  to  change  signif- 
icantly the  pressure  distribution  in  the  flow,  but  do  not  directly  influence 
the  velocity  distribution  and  the  onset  of  flow  separation.  This  is  due  to 
the  Irrotational  nature  of  the  Lorentz  force,  V x [V  x ^ x 5]  = ^,  due  to 
the  axial  (non-curved)  magnetic  field  ^(r,0). 

1 : 
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INTRODUCTION 


Analytical  similarity  solutions  for  electrically  conducting, 

incompressible  flows  between  inclined  plane  walls  ("diffuser")  exist 

in  presence  of  azimuthal  external  magnetic  fields  if  induced  magnetic 

1)  2) 

fields  are  neglected—  or  considered—  (depending  on  the  magnitude  of 

the  magnetic  Reynolds  number),  and  even  if  viscous  stress  relaxation  is 
3) 

taken  into  account—  . For  compressible  plasma  flows  in  diffusers  with 

azimuthal  magnetic  fields,  closed-form  solutions  were  obtained  in  the 

4) 

polytropic  approximation—  , whereas  the  general  case  with  heat  transfer 
and  Ohmic  and  viscous  heating  required  both  similarity  and  numerical 
methods^^^  If  the  applied  magnetic  field  b^  is  azimuthal,  the  force 
field  1 X S resulting  from  the  interaction  of  the  current  density  2 and 
magnetic  induction  b is  rotational,  V x (^  x ^)  = ^ V x (^-V^)  f 

since  ^*7^  ^ 5 (curved  ^-lines).  As  a consequence,  the  velocity  field 
V in  the  diffuser  is  modified  by  the  Lorentz  force,  whereas  the  pressure 
distribution  is  similar  to  that  of  the  (nonconducting)  Hamel  f low^^ . 
Similarity  solutions  for  magnetohydrodynamic  flows  in  diffusers  with 
homogeneous  external  magnetic  fields  have  not  been  found  previously.  We 
will  demonstrate  that  a two-dimensional  (r,0)  similarity  solution  exists 
also  for  incompressible  magnetohydrodynamic  diffuser  flows,  if  the 
applied  magnetic  field  is  homogeneous  and  in  the  axial  direction  (z). 

The  Induced  magnetic  field  is  then  also  axial  (and  z-independent)  so 
that  the  magnetic  force  density  is  Irrotational,  V x (^  x ^)  * since 
J X ^ ■ -7(^^/2p^)  + ^*7^  where  b*7b  = ^ (straight  b-llnes).  In  the 

latter  case,  therefore,  the  magnetic  volume  force  changes  the  pressure 
distribution  in  the  flow  but  does  not  affect  the  velocity  field  directly. 


4 


We  present  a similarity  transformation  which  transforms  the  magneto- 
hydrodynamic  equation  of  motion  for  the  velocity  field  v(r,0)  and  the 
Induction  equation  for  the  magnetic  field  b(r,6)  for  the  incompressible 
diffuser  flow  with  an  applied  liomogeneous  field  b^  In  the  axial  direction 


into  nonlinear  coupled  differential  equations  for  the  amplitude  fields 
f(0)  and  ii((0)  of  v(r,0)  and  b(r,0),  respectively.  The  boundary-value 
problem  for  f(0)  is  solved  in  closed  form,  whereas  the  eigen-value  prob- 
lem for  ip(0)  is  treated  by  means  of  perturbation  theory  assuming  magnetic 

Reynolds  numbers  < 1 (R^  > 1 is  generally  not  encountered  in  experi- 
8) 

ments)—  . Graphs  of  the  velocity  and  magnetic  field  are  presented.  The 
pressure  distribution  and  the  onset  of  separation  of  the  laminar  boundary 
layer  from  the  walls  in  the  magnetohydrodynamic  flow  are  discussed. 

BOUNDARY-VALUE  PROBLEM 

We  consider  the  radial  outflow  v = {u,0,0}  of  a viscous,  electrically 
conducting  incompressible  fluid  (liquids  such  as  mercury  or  subsonic, 
quasi- incompressible  plasmas)  in  a diffuser  with  plane  walls  at  0 = ±0^ 
and  a homogeneous  magnetic  field  b^  = {0,  0,  b^}  in  the  axial  direction 
(Fig.  1).  The  diffuser  walls  extend  radially  from  r^  > 0 to  r2  < “ 


and  are  ideal  conductors  which  collect  the  aximuthal  currents  which 

are  Induced  by  the  radial  flow  of  the  conducting  fluid  across  the  magnetic  | 

^ I 

field  b in  the  z-direction.  The  electrodes  are  connected  by  a load  ] 

circuit  through  which  the  current  i (amp)  is  transported.  Boundary  effectst 

at  the  remaining  insulating  walls  in  the  planes  z = ±z^  are  not  considered  | 

by  assuming  that  2z  >>  0 (r,+r„)/2,  so  that  all  flow  fields  are  two-  * 

“»  o 1 z 

dimensional  of  the  form  F » F(r,0).  Similarly,  hydrodynamic  and  electro-  ( 


magnetic  entrance  (r**rj^)  and  exit  (r“r2)  effects  are  disregarded.  In 


actual  diffuser  flows,  end-effects  are  minimized  by  operating  the  diffuser 
in  a similar,  larger  diffuser  duct  through  which  the  working  fluid  is 
pumped. 


For  fields  Independent  of  z.  Maxwell's  equation  V x b “ Indi- 
cates that  the  azimuthal  current  density  Jg(r,0)  produces  an  axial  magnetic 

field  b (r,0),  which  in  turn  requires  a radial  current  density  i (r,0) 
z r 

since 


p j.  = -9b  /9r,  p j =+r  9b  /96. 

^o-'0  z ’ Cr  z 

It  is  also  seen  that  the  fields  j = {j  ,jQ,0}  and  b = {0,0, b + b } 

r y o z 

satisfy  the  basic  equations  V*^  = 0 and  V‘b  = 0.  In  accordance  with  the 

current  transport  equation,  j = a(e  + v x b)  in  absence  of  the  Hall- 
8 ^ 

effect—  , the  electric  field  is  of  the  form  e = e^,  0}  and  satisfies 

the  basic  equation  V x e = ^.  We  Introduce  dimensionless  flow  fields  and 
independent  variables  in  accordance  with  (p  = pressure  field), 
u(p,0)  - u(r,e)/u^,  P(p,0)  = p(r,e)/p^,  jp^g(p,0)  = 


B(p,9)  = b (r,0)/b  , B = 1,  E ft(P.0)  = e (r,e)/e  , 
z oo  P,o  r,o  o 


and 


where 


“0  “ '^^’^0'°^  ^ Po  “ Po'^o’  % ° “o’^o’  Jo  ° P'^o'^o 


(1) 

. (2) 

, (3) 

. (4) 


Based  on  the  magnetohydrodynamic  equations  for  constant  density  p^, 
viscosity  p,  and  conductivity  o and  Maxwell's  equations,  we  derive  for 
the  dimensionless  fields  U(p,0),  P(p,0),  and  B(p,6)  of  the  plane  diffuser 
with  axial  external  magnetic  field  the  nonlinear  boundary-value  problem; 


6 


- _L  rp  + i Mn+R^2,  . rJ: 

^3p  3p  ^ 2 ^ ^ 2 ,,2 

p 00 


0 - [-P  - |W(l+B)^  + 2R“^ 


f H . 0 

3p  P 


3^B  , 1 3B  ^ 1 3^B  p 1 3 f 

— 7 7 77  ■*■  "T  — 7 “ ^ 7 77  [PU (I"*"®) 

3p2  P p2  30^  ^ P ^P 


with 


U(P,0) 


0=±0 


0,  Pj^  1 P 1 P2 


[3B(p,0)/30]Q^_^g  =0,  Pj^  1 P 1 P2 


, (5) 

, (f>> 

, (7) 

, (8) 

, (9) 

.(10) 


+0O 

J U(p,0)  pd0  = Q > 0 
-0 


,(11) 


/ [3B(p.0)/3p]g^^Q  dp  » -Rj^l 


.(12) 


q • p u r Q and  i = 1 r I are  the  dimensional  flow  rate  and  dimen- 
o o o o o 

slonal  load  current,  each  per  unit  length  Az  = 1.  Although  asymmetric 
Hamel  flows  are  compatible  with  the  basic  equations  and  boundary  condi- 
tlons^\  we  consider  exclusively  symmetric  outflows  with 

U(P,-9)  - U(P,+e),  B(P,-0)  = B(p,+9)  . (13) 

The  flow  Reynolds  number  R,  the  magnetic  Reynolds  number  Rj^,  and  the 
energy  ratio  N are  defined  as 

0 7 

(14) 


2 2 

000  m 000  0000 


Eqs.  (5-6),  (7),  and  (8)  are  the  equation  of  motion  (r  and  0 components) , 
the  continuity  equation,  and  the  induction  equation  (V^b  = “V^p  ^ x (^b)]. 


respectively.  The  boundary  conditions  (9)  and  (10)  consider  that  the 


i 


[ , 


fluid  does  not  slip  at  the  walls  0 = ±6  and  that  E (p,±6  ) = J (p,±0  ) “ 

o p o p o 

R ^p  ^ 3B(p,±0  )/30  “ 0 at  the  electrodes  of  infinite  conJuctivity , 
m o 

respectively.  Eqs.  (12)  and  (13)  specify  the  (given)  mass  (Q)  and 
current  (I)  flows. 

From  the  solutions  U(p,0)  and  B(p,0)  of  Eqs.  (5)-(13),  the  remaining 
(dimensionless)  flow  fields  are  obtained  as 

J = R^^  p“^  3B/30,  J.  = -R~^  3B/9p  ,(15) 

p « 0 m 

Ep  ^ 3B/30.  Eg  = U(l+B)  3B/3p  .(16) 

SIMILARITY  TRANSFORMATION  AND  EIGENVALUE  PROBLEMS 
The  boundary-value  problem  in  Eqs.  (5)- (12)  has  self similar 
solutions  of  the  form  F(p,0)  = p'^  G(0)  for  the  velocity  field  U(p,0)  and 
the  magnetic  field  B(p,6),  whereas  the  pressure  field  P(p,0)  has  a more 
complicated  structure.  Integration  of  Eq.  (6)  yields  for  the  sum  of 
hydrostatic  and  magnetic  pressures 


P + y N(1+B)^  = 2R~^  ^ + B(p) 


(17) 


where  B(p)  is  an  integration  "constant"  with  respect  to  0.  Elimination 
1 2 

of  P + -^(1+B)  from  Eq.  (5)  by  means  of  Eq.  (17)  results  in 


.(18) 


i;^  = p-lfoU  2 3U  1 3^,  dg(p) 

^ ^ - p ^ J ..2^  ■ dp 

p p 30 

1 2 

Alternatively,  V[P  + — N(l+B)  ] could  have  been  eliminated  from  Eqs. 
(5)-(6)  by  taking  the  curl  of  the  equation  of  motion.  Since  3(pU)/3p  » 0 
by  Eq.  (7),  Eq.  (8)  can  be  reduced  to 


ill  + i M + -L  iii 

3p2  P p2  302 


Rm  ^ 


i® 

3p 


.(19) 


It  is  now  seen  that  Eqs.  (18)  and  (19)  can  be  reduced  to  ordinary 


S 


differential  equations  by  means  of  the  similarity  transformation 


f ' 


► 


\ 


U(p,0)  - f(e)/p",  n - 1 
B(p,e)  “ ijj(e)/p”*,  m > 0 

where  n ••  1 by  Eq,  (7)  and  m is  a still  undetermined  eigenvalue, 
namely: 


~ + Rf^  + 4f 
de"^ 


Rp^^ 

dp 


= -aR 


I + m(m+Rj^f)i|)  = 0 


, (20) 
. (21) 


, (22) 

. (23) 


In  Eq.  (22),  oR  is  a separation  constant,  and  the  p-dependent  con- 
stituent of  the  pressure  field  becomes 

e(p)  = ^p“^  + (24) 

where  is  the  hydrostatic  overpressure.  Combining  of  Eqs.  (17)  and 
(24)  gives  for  the  total  pressure  field 


P(p,0)  + |W[l+B(p,0)]^  - 2R"^  + i -^  + P^  . (25) 

P P 

Note  that  the  sum  of  hydrostatic  and  magnetic  pressures  is  of  the  same 
selfsimilar  form  as  the  pressure  field  in  the  classical  Hamel  flow^^ . 

Eqs.  (20)-(21)  and  (25)  represent  the  similarity  transformation  for  the 
nonlinear  magnetohydrodynamic  flow  problem  under  consideration. 

Eq.  (22)  and  the  boundary  condition  (9)  give  for  the  determination 
of  the  velocity  amplitude  f(0)  the  nonlinear  eigenvalue  problem  with 
eigenvalue  a; 

.2 

— j + Rf  + 4f  + oR  - 0 , (26) 

d0 

f(±6o)  - 0 , (27) 

9 


J 


1 


with 


f(-e)  - f(+0),  f(0)  S 1 . (28) 

Eq.  (23)  and  the  boundary  conditions  (10)  and  (12)  give  for  the 
determination  of  the  magnetic  field  amplitude  i/<(0)  the  linear  eigenvalue 
problem  with  eigenvalue  y: 


+ y ij)  = f (0)<^  , y H 


d0 


di|;(±0^)/d0  = 0 


, (29) 
. (30) 


with 


Ij 

i|)(-9)  = 'j<(+e)  11/(160)  ° 


. (31) 


In  (4),  the  identity  f(0)  = 1 is  due  to  the  nondimensionali- 


zation  of  U(p,0)  with  respect  to  u = u(r  ,0),  Eq.  (4).  Since  B(p,0) 

o o 


is  non-dimens ionalized  with  respect  to  the  arbitrary  (external)  b^ < 
|i|((0)|  1 in  general  in  accordance  with  Eq.  (31).  In  the  eigenvalue 

problem  (29)-(31),  f(0)  represents  a variable  coefficient  which  is 
assumed  to  be  known  from  the  solution  of  Eqs.  (26)-(28). 

Since  iJ/(0)  is  an  even  function  of  0,  by  Eq.  (31),  dip(6)/d0  is  an 


odd  function  of  0.  Hence,  the  flux  of  the  radial  current  density  Jp(p,0) 


through  any  surface  Pj^  £ P = constant  £ p^,  |0|  0o»  vanishes  [Eq.  (15)], 


+6, 


+0 


pde  . r;,/  p-"  /_  ft-  ae  . 0 


-0 


-0 


• (32) 


SOLUTION  FOR  f(0) 

We  seek  symmetrical  solutions  of  the  boundary-value  problem  (26)-(28) 
which  represent  pure  outflows,  0 £ f(0)  £ 1.  Net  outflow  solutions  (Q > 0) 
with  backflow  regions  where  f(0)  < 0 are  not  discussed  since  they  are  in 
all  probability  not  stable  (flow  separation) . A first  integration  of 
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Eq.  (26)  gives 


I 


i 


f 


(33) 

where  [df (t9^)/d0]^  - (2/3)RC  ^ 0 for  K > 0.  Since  df(O)/d0  = 0 and 
f(0)  » 1, 

C“l+6R^+3a>0  . (34) 

Backflows  occur  at  the  walls  0 = +0  if  df/d©  < 0 for  0 = ±0  , 

o ^ o 

Accordingly,  the  eigenvalues  a of  pure  outflows  are  larger  than  a 
critical  value  which  is  negative, 

a ^a,  a = - ^(1+6R"^)  . (35) 

In  view  of  Eq.  (34),  the  trinomial  in  Eq.  (33)  has  the  roots 

^1  “ ^2  ~ ^+’  ^3  “ 

where 

f^  * ■|{-(1+6R"^)  ± [(l+6R"^)^-12a-4(l+6R"^)]^^^}  (37) 

with 


f.  . cmlagate  complex,  ^ 


1 4' 


(38) 


Formal  integration  of  Eq.  (34)  yields,  under  consideration  of  f(0) 


fl  = 1. 


(2R/3)^^^0  = ± / ' [(-I)(f-f^)(f-f2)(f-f3)]' ^'‘df 

In  the  evaluation  of  the  elliptic  integral  in  Eq.  (39)  for  pure 
outflows,  0 ^ f(0)  ^ 1,  two  cases  have  to  be  distinguished: 


f(0) 


1-1/2. 


(39) 


CASE  I:  f^  = conjugate  complex. 

In  this  case,  the  eigenvalues  a lie  in  the  region  defined  by  Eq.  (38), 

2 

The  substitution  f * X [ (l-coS(J))/ (l+cos(t))  ] , 0 ^ (|i  ^ n,  reduces 

1/2  - _i 

Eq.  (39)  to  (2R/3)  0 “ + X F(^,k).  The  solution  f(0)  is  obtained 

11 


j 


by  inversion  of  the  elliptic  integral  F(((i,k)  of  the  first  kind  as: 


f(0)  *1-  X 


2 l-cn[(2R/3)^^^Xe;k] 
l-fcn[2R/3)^'^^X9;k] 


where 


X^=  [3(1-Hj+4R 

- |[l+|x”^(l+2R"^)] 


. (41) 


. (42) 


The  eigenvalue  a is  given  by  the  boundary  condition  (27)  as  the  real 


root  of  the  transcendental  equation, 

cn[(2R/3)^^^Xe^;k]  = (X^-1)/ (X^+1) 


. (43) 


CASE  II:  f^  = real. 


In  this  case,  the  eigenvalues  lie  in  the  interval,  -(1+6R  ■*‘)/3  _< 

“ ” (1+6R  ^)  (1-2R  ^)/4,  and  ^ ^2  ^ ^3  (36)-(37)  since 

2 

f^  i O'  The  substitution  f = f2  + (f2“f2^^°®  ‘f’*  0 ^ i ''</2,  reduces 
1/2  - -1 

Eq.  (39)  to  (2R/3)  X0  = + X F((|i,k).  Inversion  of  the  elliptic  integral 
yields  the  solution: 

f(0)  = f2  + (l-f2)cn^[(2R/3)^''^X0;k]  (44) 


where 

X^  = (1-13) /4 
k^  = (l-f2)/(l-f3) 

The  eigenvalue  a is  determined  by  the  boundary  conditions  (27)  as 
the  real  root  of  the  transcendental  equation, 
cn^[(2R/3)^''^X0^;k]  - - f2/(l-f2) 


. (45) 

• (46) 


. (47) 


For  large  Reynolds  numbers,  R >>  1,  as  encountered  in  most  exper- 
iments, both  Eqs.  (40)-(43)  and  Eqs.  (44)-(47)  Indicate  that  X,  k,  and 
a - and,  hence,  the  solution  f(0/0^)  - depend  only  on  the  combination 
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1/2 

I * R 0 . This  means  that  flows  with  different  R and  0 , but  with  the 
o o 

1/2 

same  I ■ R 6 , have  Identical  velocity  distributions  f(0/6  ).  This 
o ' o 

invariance  principle  is  a special  case  of  the  one  first  discovered  for 

A) 

compressible  magnetohydrodynamlc  flows—  . 

At  the  transition  from  pure  outflow,  f(0)  ^ 0,  to  backflow  f(0)  < 0 

in  the  vicinity  of  the  walls  0 = ±6^*  the  azimuthal  velocity  gradient 

drops  to  zero  at  the  walls.  This  transition  point  is  usually  defined 

as  the  onset  of  flow  separation  in  incompressible  fluids^\  Appllca- 

of  the  condition  df(±0^)/d0  = 0 to  the  solution  in  Eq.  (44),  which  is 

valid  for  eigenvalues  down  to  the  separation  value  a = -(1+6R  ^)/3, 

gives  for  the  critical  duct  angle  0 beyond  which  separation  occurs 

o 

(2R/3)^'^^X0^  = K(k)  (48) 

9) 

where  K(k)  is  the  complete  elliptic  Integral  of  the  first  kind—  . 
Substitution  of  Eq.  (48)  into  Eq.  (47)  shows  that  f2  = 0 so  that  at  the 
onset  of  separation 

fj^  = 1,  f2  = 0,  f^  - -(1+6R"^)  (49) 

and 

= (l+3R“^)/2,  k^  = (1+3R"^)"^/2  (50) 

s 

Combining  of  Eqs.  (48) -(50)  shows  that  the  critical  duct  angle  0^ 
beyond  which  (0^  > 0^)  flow  separation  sets  in  decreases  with  increasing 
Reynolds  number  R, 

Qq  “ [3/(3+R)r^  K(^'R/(6+2R))  . (51) 

Since  the  eigenvalue  a «■  -(1+6R  ^)/3  is  negative  at  the  onset  of 

1 2 

separation,  the  total  wall  pressure  P(p,±0^)  + •^[1+B(p ,±0^) ] is 
positive  at  the  onset  of  separation  by  Eq.  (25)  if  the  overpressure 
is  sufficiently  positive. 

In  order  to  have  a stable,  unseparated  diffuser  flow,  experiments 


are  generally  operated  at  sufficiently  positive  eigenvalues,  a ^ -(1+6R  )/3. 

Fig.  2 gives  the  eigenvalue  a for  pure  outflows  versus  the  invariance 
parameter  I ■ 0^  (R  >>  1)  based  on  Eqs.  (43)  and  (47)  for  the  solu- 

tions of  types  I and  II,  respectively,  a decreases  monotonically  from 
large  positive  values  (for  large  1-values)  to  its  minimum  value  a = - 1/3 
(for  I “ 3.18704)  below  which  backflows  occur  at  the  walls. 

Fig.  3 shows  the  azimuthal  distribution  f(x),  x = of  the  radial 

flow  velocity  of  pure  outflows  for  eigenvalues  a = - 1/3  to  5.0  (I  = 3.187 


to  0.591)  based  on  Eqs.  (40)  and  (44)  for  the  solutions  of  types  I and  II, 
respectively.  It  is  seen  how  f(x)  changes  from  the  shape  of  regular  out- 
flows for  small  values  of  I = 0 to  the  critical  flow  with  the  maximum 

o 

value  I = 0^  = 3.187  (of  a pure  outflow),  at  which  the  transition  to 

mixed  flows  occurs.  The  curves  in  Fig.  3 imply  large  Reynolds  numbers, 

R >>  1,  for  which  the  Invariance  principle  holds. 

Fig.  4 shows  the  critical  duct  angle  0^,  beyond  which  (6^  > 0^) 
backflows  occur,  versus  the  Reynolds  number  R based  on  Eq.  (51),  It  is 
recognized  that  0^  decreases  rapidly  with  increasing  R.  Pure  (unseparated) 
outflows  can,  therefore,  only  be  achieved  at  large  Reynolds  numbers  R 
if  the  duct  angle  is  sufficiently  small,  0^  ^ 0^(R).  It  should  be  noted 
that  in  the  present  case  of  a homogeneous,  axial  b^,  the  external  mag- 
netic field  can  not  prevent  flow  separation,  in  contrast  to  an  azimuthal, 

external  magnetic  field  ^ (r)  which  inhibits  flow  separation  for  Hartmann 

o 

2 A ^ 

numbers  H ^ 2 R/3-—  . The  physical  reason  for  this  is  to  be  seen  in 

the  different  nature  of  the  Lorentz  force,  which  is  irrotational 
[7  X (T  X b)  - in  the  former  case  and  rotational  [V  x (j  x b)  in 

the  latter  case. 
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SOLUTION  FOR  Xe) 

For  magnetohydrodynamlc  flows,  the  magnetic  Reynolds  number  is 
commonly  not  large,  - 1.  For  this  reason,  the  eigenvalue  problem 
(29)-(3l)  for  the  amplitude  ij<(6)  of  the  Induced  magnetic  field  B(p,0) 
can  be  solved  by  means  of  perturbation  theory,  assuming  that  is 

a small  term  in  Eq.  (29).  Since  f(6)  ^ f(0)  = 1,  this  assumption  is 
valid  for 


« ■7‘’> 

n 


n = 1,2,3,. . . 


In  the  O-th  approximation,  ->  0,  the  eigenvalues  and  eigenfunctions  of 
Eqs.  (29)-(30)  are: 


(nir/e^)  , n = 1,2,3,. . 


, (53) 


■'(9)  = c cos(mT9/0  ),  n = 1,2,3,... 
n n o 


. (54) 


with  •^(0)  “ 0 for  n = 0 by  Eq.  (31).  For  the  determination  of  the 
eigenvalues  and  eigenfunctions  of  the  perturbed  problem  (29)-(31), 

Rj^  > 0,  we  set 

, (0)  (1) 

00 

’I'  (e)  = i^^°^(0)  + I 1)'^°^(0)  + ... 

n ’^n  , vn  V 

v=l 

in  accordance  with  1-st  order  perturbation  theory.  Substitution  of 
Eqs.  (55)-(56)  into  Eq.  (29)  yields  in  the  same  approximation 


VW  ' .<»  v"»  ♦<»>  + ... 

n M ’^n  , vn  V v 


^n  ^n  ^n  , vn  ^v 
v«l 

since  the  terms  of  0-th  order  drop  out.  Multiplication  of  Eq.  (57)  by 


15 


1 


(0) 


and  Integration  over  (“9^^.  + 6^)  yields  for  the  perturbation  of 


the  eigenvalues 


+0 


(1) 


ynn/6  ) f cos  (mi0/6  )f(0)d(0/0  ) 
o'.  o o 

— 0 

o 


, (58) 


upon  consideration  of  Eqs.  (53)-(54).  Multiplication  of  Eq.  (57)  by 


(0) 


, V n,  and  integration  over  (-9^.  + 9^)  yields  for  the  expansion 


coefficients  of  l-st  order 


+0 


(1)  ^ ° 

= 9 9 / cos(vTr0/0  ) f (0)cos(mT0/0  )d(0/0  ) 


vn 


(vTr/0  ) -(mr^  ) v -0 
o o o 


. (59) 

upon  consideration  of  Eqs.  (53)-(54).  Combining  of  Eqs.  (58)-(59)  and 
(55)-(56)  yields  for  the  eigenvalues  and  eigenfunctions  of  Eqs.  (29)-(31): 

+1 


^n  “ (nir/e^)  [ (mT/0^)  - / cos  nirx  f(x)  dx],  n = 1,2,3,.. 


(60) 


- ^ 

-’‘2  > 


[cos(mT0/0  )+R^(mT/0  ) Z'  b cos(mr0/0  )] 
oM  o.vn  o 

v=l 


[(-l)"+R^(mi/0^)  Z'  (-I)''  b^^] 


, n=l,2,3,., 


(61) 


v=l 


where 


+1 


f 2 2 

b^^  = j cosvnx  f(x)  cosniTx  dx/[(viT/0^)  -(mr/0^)  ] 


-1 


(62) 


In  Eqs.  (60)-(62),  x = previously,  and  the  expansion  coeffi- 


cients c have  been  determined  in  Eq.  (61)  by  means  of  Eq.  (31).  The 


parameter  m,  which  determines  the  radial  decay  of  the  Induced  magnetic 
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. (63) 


field  B(p,6)  • p '|'(0),  is  obtained  from  p , 

n 


m = m 


Since  an  infinite  number  "n"  of  eigenfunctions  '<'^^(0)  exists 
mathematically,  the  solution  for  the  induced  magnetic  field,  B(p,0)  = 
p is  not  unique.  Probably,  only  the  eigenfunction  with 

n“l  is  physically  stable,  and  can  be  observed  in  actual  diffuser  flows 
across  an  axial,  homogeneous,  external  magnetic  field.  A similar  non- 
uniqueness exists  for  the  velocity  field  U(p,0)  = f(0)/p  of  the  Hamel 
flo»^\  For  one  and  the  same  flow  rate  Q > 0 and  duct  angle  0^,  up  to 
an  infinite  number  of  eigenfunctions  f(0)  exist  mathematically,  which 
represent  pure  outflow  (f>0),  symmetrical  mixed  out-  and  inflows  (f^O), 
and  asymmetrical  flows  (f^O)— ^ . The  solution  for  f(0)  established  above 
is  only  insofar  "unique"  as  we  have  excluded  all  types  of  mixed  flows 
(f<0)  for  given  values  of  Q > 0 and  0^  by  means  of  physical  arguments. 

In  Table  I,  the  first  few  (n*l,2,3)  eigenvalues  m^  are  compared 

with  their  values  m^^^  in  the  0-th  approximation  for  0 = 5°,  a = -1/3, 

n o 

5,  and  = 10  10  10^  based  on  Eqs.  (60)  and  (63).  It  is  obvious 

that  m differs  hardly  from  m^^^  and  varies  very  little  with  a and  R,^, 
n n m 

as  long  as  R^^  <<  rnr/d^. 


TABLE  I:  Eigenvalues  " ^n  ~ 1»2,3(0^=5°) . 


R=10"^ 

a— 1/3 

R»10"^ 

R-10° 

R-IO"^ 

0*5 

R-IO"^ 

R=10° 

n 

(0) 

”l 

35.9976 

35.9760 

35.7589 

35.9969 

35.9694 

35.6927 

1 

: 36.0000 

“2 

71.9977 

71.9770 

71.7695 

71.9968 

71.9675 

71.6748 

2 

: 72.0000 

"3 

107.9977 

107.9770 

107.7696 

107.9967 

107.9672 

107.6717 

3 

: 108.0000 
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Flf^s.  5 and  6 show  the  first  few  (n=l,2,3)  eigenfunctions  4;^(x)/Rj^I 

versus  x ■»  9/9^  for  a *•  -1/3,  = 10°  and  a = 5,  = 10  respectively. 

In  both  cases,  it  is  assumed  that  0 =5°  and  p,  = 1 (r,  = r ) and 

o 1 I o 

P2  “ 100(r2  “ lOOr^) . It  is  recognized  that  '('jj(x)  varies  little  with  a am 
R^,  and  that  , 

Thus,  a homogeneous,  external  magnetic  field  in  the  axial  direc- 
tion has  the  following  effects  on  the  radial  flow  in  a diffuser  with 
electrodes  in  the  planes  9 = ±9^.  The  radial  motion  of  the  conducting 
fluid  across  b^  induces  radial,  J^Cp.B),  and  azimuthal,  Jg(p,0),  current 
densities.  While  the  net  radial  current  flow  due  to  Jp(p,0)  is  zero 

[Eq.  (32)],  the  net  azimuthal  current  flow  due  to  J (p,0)  is  I 0 

y 

[Eq.  (12)1  if  the  external  circuit  is  closed.  Jg(p,0)  = 
changes  its  direction  at  those  coordinate  values  0 where  Jp(p,0)  = 

1^1  p (nri-1)  assumes  an  extremum  [Eqs.  (15),  (21),  (61)]  as 

illustrated  in  Figs.  5-6.  The  external,  b^,  and  induced,  ?(p,0), 
magnetic  fields  in  the  axial  direction  produce  an  irrotational  Lorentz 
force,  which  changes  the  pressure  distribution  in  the  flow  [Eq.  (25)], 
but  has  no  direct  influence  on  the  onset  of  flow  separation  [Eq.  (51)]. 

The  selfsimilar  solutions  derived  are  based  on  the  assumption  that 
the  flow  is  two-dimensional  (r,0),  which  requires  a quasi-infinite 
extension  of  the  flow  in  the  z-direction.  Whether  the  selfsimilar  solu- 
tions are  realistic  enough  remains  to  be  shown  by  comparison  with  actual 
magnetohydrodynamic  diffuser  flow  experiments,  since  the  condition  of 
similarity  may  impose  restrictions  on  the  generality  of  the  solutions 
obtained. 

The  radial,  two-dimensional  diffuser  flow  considered  represents 
an  idealization  which  is  analogous  to  the  fully  developed,  axial 


1 


) 
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Hagen-Polseuille  flow  between  parallel  walls.  For  a more  realistic 
theory,  hydrodynamic  and  electromagnetic  end-effects  would  have  to  be 
considered,  as  well  as  z-dependence  of  the  flow  fields  due  to  insulating 
walls  in  the  planes  z “ ±z^,  < <».  These  complications  would  render 

it  extremely  difficult  to  find  self similar  solutions,  presumed  that 
similarity  transformations  exist  in  the  three-dimensional  case  at  all. 
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Similarity  analysis  of  magnetohydrodynamic  flows  with  viscous  stress  relaxation 

H.  E.  Wilhelm  and  S.  H.  Choi 

Dtpartmvnt  of  Electrual  Engineering.  Colorath  State  University,  i'ort  Collins.  Colorado  H0S2.t 
(Reccivctl  I March  1977) 

A novel  similarity  solution  in  terms  of  a hyperelliptic  integral  is  given  for  a magnetohydrodynamic  flow 
across  an  azimuthal  magnetic  field  in  a diverging  duct,  under  consideration  of  viscous  stress  relaxation 
Velocity  profiles  and  the  critical  duct  angle  for  flow  separation  are  calculated  as  a function  of  the  Reynolds 
number  and  the  Hartmann  number.  It  is  shown  that  viscous  stress  relaxation  modifies  the  velocity 
distribution  and  reduces  considerably  the  critical  duct  angle  for  flow  separation  at  low  Reynolds  numbers  At 
large  Reynolds  numbers,  viscous  stress  relaxation  is  less  important,  and  the  results  approach  asymptotically 
those  of  ordinary  magnetofluiddynamics,  which  is  based  on  a static  relation  between  viscous  stresses  and  the 
velocity  component  gradients. 


I INTRODUCTION 

In  classical  fluid  mechanics'  and  magnetohydro- 
dynamics,‘ it  is  assumed  that  Inhomogeneities 

in  the  velocity  components  I’y  produce  instant- 
aneously viscous  stresses  n,y.  Mathematically, 
this  is  expressed  through  a phenomenological 
‘'flux”-"force”  relation,  given  for  incompressible 
fluids  or  subsonic  flows  by*'* 

In  a real  continuum,  velocity  inhomogeneities  do 
not  switch  on  viscous  stresses  instantaneously  but 
rather  in  accordance  with  a relaxation  process 
of  characteristic  time  t.  indeed,  the  moment 

II, (c,c,-5e’6„y(c,r,f)d'c 

of  the  Boltzmann  equation  yields,  for  an  incom- 
pressible viscous  fluid  (V»v  = 0),  the  relaxation 
equation  (discussed  in  connection  with  the  13-mo- 
ment approximation  in  the  Appendix); 

®»n,y+  -’■“n,,  -/)(V,y^+ 

if  thermal  forces  and  terms  of  higher  order  in  the 
field  derivatives  are  neglected.  The  viscosity  p, 
the  fluid  pressure  p,  and  the  viscous  stress  re- 
laxation time  T are  interrelated  by  tJ.=pr.  This 
equation  satisfies  the  basic  requirements  of  a 
(classical)  r, /-dependent  field  equation  since  it 
(i)  contains  space  and  time  derivatives,  and  (ii) 
is  invariant  against  Galilei  transformation.  The 
static  Navier-Stokes  relation  is  not  in  accordance 
with  either  requirement. 

The  third  equation  is  in  the  typical  form  of  an 
Inhomogeneous  relaxation  equation,  with  a forcing 
term,  -p{'^,Vj+  ^/t’,).  It  is  seen  that  it  reduces 
to  the  first  equation  if  the  temporal  (8,n,y)  and 
convective  (t^,^,n,y)  relaxation  terms  are  dis- 
regarded. 


According  to  the  first  equation,  viscous  stresses 
would  propagate  in  accordance  with  a (parabolic) 
diffusion  equation  (continuous  "signals”  and  in- 
finite speed  of  propagation).  According  to  the 
third  equation,  viscous  stresses  would  propagate 
in  accordance  with  a (hyperbolic)  wave  equation 
(discontinuous  "signal”  and  finite  speed  of  propa- 
gation). This  is  readily  shown,  e.g.,  by  combining 
the  equation  of  motion  for  the  viscous  fluid  with 
the  first  and  third  stress  transport  equations, 
respectively,  for  the  case  of  a small  one -dimen- 
sional velocity  perturbation.  Thus,  these  equa- 
tions give  rise  to  a qualitatively  significant  dif- 
ferent behavior  in  viscous  stress  transport. 
Quantitatively,  the  term  9,n,,  is  of  importance 
for  short  processes  with  a duration  time  / -s  t = p / 
p.  A criterion  for  the  quantitative  significance  of 
the  term  is  not  as  easily  establishable , 

since  ^,n,,  may  be  quasisingular  at  certain  points 
of  the  fluid.  Similarly,  the  rigorous  theory  of  heat 
transport  has  to  be  based  on  a (hyperbolic)  wave 
equation.' 

We  consider  herein  subsonic  flows  of  dense, 
ionized  gases  across  an  external  azimuthal  mag- 
netic field  in  a duct  with  inclined  walls  (so- 
called  diffuser.  Fig.  1).  The  analysis  is  based  on 
the  magnetohydrodynamic  equations  with  viscous- 
stress  relaxation,  i.e.,  we  disregard  effects  of 
"magnetic”  viscosity  (which  occur  in  highly  rare- 
fied plasma  flows)  assuming  that  u>,r,  « 1,  where 
w,  = e,  Bg/m,  and  t,  are  the  gyration  frequency  and 
collision  time  of  the  ions,  respectively.  By  means 
of  an  exact  (nonlinear)  similarity  solution,  we 
demonstrate  that  convective -stress  relaxation 
affects  the  onset  of  flow  separation,  i.e.,  the 
first  occurrence  of  wall  back  flows  which,  in  gen- 
eral, are  unstable  and  result  in  a turbulent  bound- 
ary layer.  Flow  separation  is  commonly  observed 
if  for  given  Reynolds  (A)  and  Hartmann  (3C)  num- 
bers, the  duct  angle  6g  is  increased  beyond  a 
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critical  value  0^.  The  calculated  velocity  distri- 
butions are  qualitatively  in  agreement  with  velo- 
city profiles  observed  in  diffusers/ 

Magnetohydrodynamic  diffusers  with  transverse 
magnetic  fields  and  nonvanishing  electric  load  are 
frequently  used  to  study  the  transformation  of 
kinetic  flow  energy  (due  to  thermal  expansion)  into 
electric  energy.  For  the  experimental  realization, 
it  is  suitable  to  install  the  diffuser  in  a similar 
larger  diverging  duct  through  which  the  working 
fluid  is  pumped  at  a constant  rate  in  order  to  mini- 
mize three-dimensional  entrance  effects.  The 
development  of  the  boundary-layer  with  entrance 
effect  is  a complex  problem  which  has  been  ana- 
lyzed only  for  magnetohydrodynamic  flows  between 
nonlncllned  walls  by  means  of  Goertler  series  ex- 
pansions.’ 

II  NONLINEAR  BOUNDARY-VALUE  PROBLEM 

Let  cylindrical  coordinates  (r,  0,z)  be  introduced 
for  the  description  of  the  magnetohydrodynamic 
flow  model  (Fig.  1).  The  conducting  fluid  is  bound- 
ed in  the  surfaces  (9=-f  B^.rjSrSrj)  and  (e=~9^, 
r,  ^r^r,)  by  insulating  walls,  and  in  the  surfaces 
(z  B-f  2 .)  and  (2  B -2.)  by  electrodes,  which  are 
connected  through  an  ideal  circuit  (R  ■ 0).  The 
conducting  fluid  is  injected  through  the  inner  cyl- 
inder section  (rBr,,- 

and  removed  downstream  through  the  outer  cylin- 
der section  (rBr„  - Ss+  63,-2  .s  9^+*.). 
The  boundary  layers  at  the  electrodes  are  disre- 
garded compared  with  those  at  the  insulating  walls 
by  assuming  that  the  interelectrode  spacing  is 
large,  + '^he  magnetic  field  has 

its  sources  in  an  electric  current  I flowing  through 
a conducting  rod  (OSrSrg,  -•>£z^+“>,rg<r^). 

In  accordance  with  Stokes’s  law,  ^B*<fsB  (igt,  the 
magnetic  field  is  azimuthal  (Uo  the  permeability 
of  vacuum)  and  has  the  induction 


that  the  Hall  effect  is  negligible  (i*»,T,«  1), 

)’^a(£,+  uB,)e,. 

The  resulting  Lorentz  force  density  is  a purely 
radial  field  which  opposes  the  inducting  flow, 

j X Bb  -(t(£,+ 

Because_of  V x E = 0 and  V • f = ^(V  • E + B • V x v 
-v*Vxb)bctV»EbO,  the  axial  electric  field  is 
inhomogeneous,  and  vanishes, 

E 0 , RbO, 

since  the  load  of  the  external  circuit  is  zero.  In 
this  case,  the  current  in  the  external  circuit  as- 
sumes the  maximum  value 

r'a  r**o 

J*o  uB^drdS. 

■'’■i  -'-•0 

These  equations  are  based  on  the  assumption  that 
the  induced  magnetic  field  is  small  compared  with 
the  external  magnetic  field,  which  implies  small 
magnetic  Reynolds  numbers,’ 

Rb‘  MoOw(»’.0)r  « 1. 

In  this  elementary  radial-flow  model,  fluid  dynam- 
ic and  electric  end  effects  at  vBr,  and  »'  = rj(|  0\ 
s 63)  are  disregarded. 

The  magnetohydrodynamic  diffuser  flow  under 
consideration  is  described  by  the  nonlinear  bound- 
ary-value problem  for  the  radial  velocity  [u 
= «(>■,  ®)],  stress  [n,yB  n,y(r,  6)],  and  pressure 
[/>■/>(»'.  ®)]  fields: 


B*  (fig/2ir)(l/r}e,,  rgSrK’^. 

The  radial  flow  VBue,  of  the  conducting  fluid 
across  the  magnetic  field  B induces  axial  electric 
(£,)  and  current  density  (j,)  fields,  presuming 


A 


ETC.  1.  Geometry  of  diverging  duct  with  radial  velo- 
city field  7 and  azimuthal  magnetic  field  Sq. 


where 


(3) 


u 8n„  _ « 

(4) 

u 8n^  „ 1 8tt 

M r-®  + n«,B  _n  _ _ 

p 8r  " r B$  ’ 

(5) 

u 8n„  _ - M 

M r + n„B  _2fj  - , 

p 8r  r 

(6) 

and 

u{r,0-±9g)~O, 

(7) 

pi  u(r,0)rd9*Q. 

(8) 

'0 
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Equation  (8)  specifies  the  flow  rate  y through  the 
diverging  duct.  For  similarity  reasons,  Eq.  (8) 
is  equivalent  to  inlet  (r  = r,)  and  outlet  (r^r,) 
boundary  conditions.  Instead  of  Eq.  (8),  it  Is 
more  convenient  to  assume  the  Reynolds  number 
R(0)  of  the  central  streamline  to  be  given," 

9»0)r/M«R(0).  (9) 

Equations  (1),(2)  are  the  r and  9 components 
of  the  equation  of  motion  of  the  conducting  fluid 
in  the  azimuthal  magnetic  field  [B,=  B,(r„)r5'Vl. 
Equation  (3)  represents  the  continuity  equation 
for  Incompressible  radial  flow,  and  Eqs.  (4)-(6) 
describe  the  convective  stress  relaxation  with  a 
viscosity  p»/>t.  It  is  noted  that  for  pure  radial 
flow 

v*vn  = «an/8r. 


III.  SIMILARITY  TRANSFORMATION 
The  dimensions  of  the  flow  fields  and  the  fluid 
constants  are  interrelated  by 

(«)“(m/pk),  (/>)•  (n,y)«  (mVpt*). 

Accordingly,  we  try  to  reduce  the  partial  Eqs. 
(l)-(6)  into  ordinary  differential  equations  by 
means  of  the  similarity  transformation; 


u{r,  9)=(M/p)r'V{9), 

(10) 

p(r,e)=(2MVp>r-»P(e), 

(11) 

n„(r,9)=(2MVp)r-*^'„(e), 

(12) 

n^{»-,9)-(2MVpVW«). 

(13) 

n„(r,  9)=  (2MVpV*VM(fl). 

(14) 

The  functions /(8),  P{8),  andg,^(8)  are  nondimen- 
sional.  Upon  substitution  of  Eqs.  (10)-(14),  we  ob- 
tain from  Eqs.  (l)-(9)  the  ordinary  boundary- 
value  problem; 


/»--4P+2g;,+  3CV, 

(15) 

P-^Pe-gH^ 

(16) 

[(P-f)/P]8rr‘f, 

(17) 

[(P-f)/P\gH‘-lf, 

(18) 

((P-/)/P] 

(19) 

where 

(20) 
(21) 

(22) 
(23) 


/(9-±9o)-0. 

r**o 

f(S)dB~Q/n, 

-^-*0 

I 

P„-P(e-±9„),  /(O)-B(O), 


From  the  stress- relaxation  Eqs.  (17)-{18)  one  ob- 
tains the  conventional  static  stress  relations  for 
^ 1/ 1 • P»  j./  j . 

IV  < lUSI  1)1  OKM  SOLU1ION 
Substitution  of  F.q.  (lO)  into  Kq.  (19)  gives  the 
pressure  function  /'(9)  in  terms  of  the  velocity 
function  /(9), 

(24) 

The  minus  sign  of  the  square  root  is  not  applicable 
as  one  verifies  by  means  of  Eqs.  (20)  and  (22).  By 
eliminating  P{9)  by  means  of  Eq.  (16)  and  g'g,  by 
means  of  Eq.  (18),  we  find  from  Eqs.  (15)  and 
(20)  the  nonlinear  boundary -value  problem  for  the 
velocity  function  f(9), 

/12/+P„+  (4/*+ />*)*/=*]+  2P„(4/»+i>?)-‘ V’* 

+ lPo-<-  (if^+Piy  (4  -3C*)/ 

+ 2[P<,+  (4/»+P^)^'*]}=0,  (25) 


where 


/(9=±9„)=0.  (26) 

The  differential  equation  for  the  corresponding 
flow  without  stress  relaxation’  is  obtained  from 
Eq.  (15)  as 

/''+/2+ (4 -3C’y+ 4P„=  0 for  P » 1/ 1 . 


The  substitution, 

^ 1^ 


(27) 


transforms  Eq.  (25)  into  a nonlinear  differential 
equation  of  first  order  for  ip- Hf), 


■^</-*+P(/)i/»+C(/)=0. 


(28) 


where 

Fi/)^  4i’o{l2/+  ^>0+  (4/’+  P?)'  '’)(V’+  PlY 

>0,  (29) 

G(/)s2[P„+{4/*+P’)‘'’] 

X {/ » + (4  - .lC*)f  + 2[P„  + ( V ’ + i*?)'  '*]} 
x[2/+P„+(4/»+P»)‘'»J-> 

SO.  (30) 


The  general  solution  of  Eq.  (28)  is  found  by  the 
method  of  variation  of  the  integration  constant  of 
the  solution  of  the  associated  homogeneous  equa- 
tion as 


0 


... — 
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4’H/)’exp(^-  F(f)d^ 

F(f)d^G{f)df\,  (31) 

/o-/(0), 

where 

-0  (32) 

for  symmetrical  flows.  Combining  of  Eqs.  (27) 
and  (31)  yields  an  integral  solution  for  9"  9(f) 
from  which  one  obtains  the  analytical  solution 
/= /(e)  by  Inversion 


rfit) 

±e-  / H-Hf)d/,  fo^m,  (33) 

Vo 

where 

«(/)-[ exp(-  ('F(f)df) 

x[^>-  jT'expJ^  jfV(/)d/)G(/)d/]}‘'*. 

(34) 

Equations  (33)  and  (34)  represent  a closed-form 
solution  for  the  magnetohydrodynamic  diffuser 
flow  with  viscous  stress  relaxation  in  terms  of  a 
hyperelllptic  integral,  since  by  Eq.  (34), 


j 


//»(/)“ 


l2/+  (V’  + P?^ 


-5(!(/^  -/?)+  5(16  _ 2K^  -P„)(/»  -/*o)  + i',;iC*(/-/„) 


+ M(4/’+-PS)’'*-(4/5+PS)’'*l 

+ 5(8  -3(?)[/(4/*  + pp*'»  -/o(4/o  + i’o)‘''] 

+ i(8  -3C»)P>{ln[2/+  (4/»+P5)‘'*]  -ln[2/o+  (4/*+ Pj)* '*]})). 


(35) 


( 


-i 

j 


I 

1 

i 

i 

i 


In  Eq.  (33),  the  ± sign  has  to  be  used  depending 
on  whether  d/ZdesO,  or  eso  In  case  of  pure  out- 
flows, f(9)>Q.  The  integration  constant  Is 
determined  in  Eq.  (32)  for  symmetrical  flows 
with  an  extremum  at  9=  0,  which  are  of  main 
practical  interest.  The  remaining  integration 
constant  contained  in  the  solution  of  Eq.  (33) 
is  determined  by  the  boundary  condition  in  Eq. 
(26),  which  gives 

/'°/f->(/)d/.  (36) 

Based  on  Eqs.  (33)-(36),  velocity  distributions 
/(S)»0  of  pure  outflows  have  been  computed  for 
the  typical  duct  angle  9j=  5°  and  given  Reynolds 
numbers  R aR(0)»/(0)  of  the  central  stream  line 
9k  0,  with  the  Hartmann  number  3C  as  parameter, 
,TC*  2 0.7R  > |R  . In  the  presence  of  viscous  stress 
relaxation,  the  onset  of  flow  separation,  as  will 
be  shown,  is  Inhibited  at  large  Reynolds  numbers 
R for  Hartmann  numbers 

» fR  for  R » 1. 

The  velocity  distributions  in  Figs.  2,  3,  and  4 
represent  net  outflows  without  backflow  regions 
since  3C>3C„,.  Figure  2 shows /(9)  for  the  rela- 
tively small  Reynolds  number  R = 10’and3C*  = 0.7R 
- 10  R(9k  5“).  The  velocity  distributions  become 
flatter  and  the  velocity  gradients  at  the  walls 
9k±  9g increase  in  magnitude  as  3C  increases. 


Figure  3 shows  /(9)  for  the  moderate  Reynolds 
number  R <=  10^  and  = 0.7R  - 5R(9„=  5°).  In  this 
case,  \df(9K  ±9„)/de\  decreases  with  decreasing 
3C  so  that  a well  developed  flow  exists  only  in  the 
central  region  for  small  3C’>  |R.  Figure  4 shows 
f(9)  for  the  relatively  large  Reynolds  number 
R = 10*  and  3C*  = 0.7R  - 5R  (9^=  5°).  For  smaU  Hart- 
mann numbers,  3C*>fR,  the  flow  is  considerably 
depressed  in  the  extended  regions  adjacent  to  the 
walls  so  that  Reynolds  numbers  R(9)a  lO*  are 
realized  only  in  the  limited  central  section  | 9| 

9 of  the  duct.  The  f(6)  curves  in  Fig.  4 show 
clearly  the  transition  to  the  limiting  velocity  dis- 
tribution, for  which  |d/(9« ±9o)/</9|  assumes  the 
smallest  realizable  value , as  3C  - |r  . It  is  con- 


FIG.  2.  /(9)  vs  9 tor  R = lO’  and  various  X'  (90=  5'). 
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eluded  that  well  developed  velocity  distributions 
exlat  for  su/flclently  large  Hartmann  numbers 
TC-TCyi;  0„),  e.g,,  3C>>  U for  ft  - 10’  (Fig.  2)  and 
TC«  ^ 2/1  lor  fl  - 10’  (Fig.  3)  and  « - 10’  (Fig,  4). 

Velocity  distributions  of  the  typical  form  shown 
In  Figs.  3 and  4 have  been  observed  In  magneto- 
hydrodynamic  diffuser-flow  experiments.’  A 
qualitative  comparison  Is,  however,  not  possible, 
since  a conical  diffuser  was  used  in  the  experi- 
ments.’ Velocity  distributions  without  stress 
relaxation  have  been  calculated  previously.' 

V.  MACNETOHYDROOYNAMIC  FLOW  SEPARATION 

The  integral  solution  derived  in  Eqs.  (33),  (35), 
and  (30)  describes  physical  flows  as  long  as  P„ 

• P{0*  ±9^)^0  lor  the  given  parameters/?,  K,  and 
0g,  since  the  pressure  field  has  to  be  positive 
everywhere,  P{0)^O,  | e|  s 0^,  The  boundary- 
value  problem  In  Eqs.  (25)-(26)  becomes.  In  the 
limiting  case  Pg^O, 

/'+J(8-3C'y  + i/?/»-0,  (37) 

/(«-*«,)- 0,  (38) 

where 

/(»-0)*l;  /?»/?(0).  (39) 

Since  0 in  Eq.  (37),  it  cannot  be  reduced  to  the 
relaxation -free  case  {P,Pg»/).  The  solution  for 
the  limiting  flow  with  vanishing  wall  pressure  Pg 
Is  by  Eqs,  (37)-(38), 

J1  <«> 

where 

«(/)-/’+ ^(8 -3C')/«-.^  (8 -3C)-1  (41) 

is  a trinomial  in  / which  has  one  real  and  two 
complex  conjugate  roots  since.  In  general. 


2/i>8  -3C’: 


A-i. 


(42) 

Equation  (40)  contains  an  elliptic  integral  which  is 
resolved  by  means  of  the  substitution 


which  gives 

±9- -(3/R)‘'*X*‘P(<(.,fe)  (44) 

where 

X « {3[  1 + (8 -.1C')«  *•]}>'’,  (45) 

*'“Hl+l^'’(2+(8-K')ft-‘]}.  (46) 


Inversion  of  Eq.  (44)  and  substitution  of  0 in  Eq. 
(43)  yields  the  explicit  solution  for  the  case  Pg 
-0; 


/(9)-l-X» 


1 -cn|[(iP)‘^X9,fe] 
l + cn|(iR)>'»xe,*]  • 


(47) 


According  to  Eq.  (44),  the  critical  duct  angle 
0g(Pg*O)  at  which  the  wall  pressure  is  zero  is 
given  by  the  boundary  condition /(9=  ±0g)^O  as, 

9„(P„.0).(3/R)‘'»X-‘P(<>)„,fe),  (48) 

d>o*arccosHx’ - ])/(X'+ l)J,  (49) 

in  terms  of  the  characteristic  flow  numbers  R and 
X [X-X(A,3C),/;=fe(P,3C)  by  Eqs.  (45)-(46)].  Since 
/^3«  0 for  3C*  - 8«  iR  by  Eq.  (42),  the  integral  in 
Eq.  (40)  diverges  for  7(9=  ±9o)=0.  Similarly,  Eq. 
(47)  diverges  in  this  case  since  fe-  1 for  3C  - 8 
- \R  by  Eq.  (46).  It  is  recognized  that 

Os9oS9o(Po«0)«irfor  3C*a8+f/?,  (50) 

0<9o£9o(Po-0)<ir  for3C’<8+|A.  (51) 


Accordingly,  physical  flow  solutions  with  P(9) 
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2 0,  I 9|  s 0^^  exist  for  all  duct  angles  0<  Og£v  if 
the  Hartmann  number  is 

KU  = B+‘,/i.  (52) 

On  the  other  hand,  if  X'-  X’^,  physical  flow  solu- 
tions with  ffO)  • 0, 1 9|  • 9g,  exist  only  for  duct 
angles  O,,--  0„(/'„« 

In  Fig.  5,  the  critical  duct  angle  fl„’  9„(/''„«0) 
for  vanishing  wall  pressure  is  plotted  versus 
H /?(0)  with  X' as  a parameter.  It  is  seen  that 
6‘  decreases  with  increasing  R,  but  increases 
with  increasing  W.  The  stabilizing  effect  of  the 
magnetic  field  at  sufficiently  large  Hartmann  num- 
bers 1C  is  apparent,  in  particular,  in  the  regions 
X^XJR). 

The  critical  duct  angle  eu(/'o=0)  is  also  ob- 
tained from  the  condition  df{8=  ±9g)/ilS^  0,  if  vis- 
cous-stress relaxation  is  not  taken  into  consider- 
ation.’ The  corresponding  curves  = Sg{R,X)  are 
shown  dashed  in  Fig.  5.  Comparison  indicates 
that  for  the  same  R and  K,  9’  "without  stress  re- 
laxation" is  considerably  larger  than  9‘  "with 
stress  relaxation”  for  relatively  small  Reynolds 
numbers,  R<10*.  Accordingly,  viscous-stress 
relaxation  has  a destabilizing  effect  on  the  flow, 
which,  however,  is  completely  negligible  tor 
large  Reynolds  numbers,  R "•  10*.  As  the  wall 
pressure  drops  to  zero,  the  laminar  flow  solution 
can  no  longer  be  realized,  and  flow  separation 
sets  in  for  duct  angles  9„>  dg(R,K)<x , 

In  the  classical  similarity  theory  for  incompres- 
sible viscous  flow  between  inclined  walls,’  solu- 
tions with  a (positive)  homogeneous  overpressure 
Po  exist  so  that  one  does  not  have  to  be  concerned 
about  negative  pressures  in  back-flow  regions  at 
the  onset  of  separation.  The  similarity  analysis 
of  the  corresponding  compressible  flow’  no  longer 
permits  solutions  with  over  pressure,  and  the 
limiting  flow  solution  with  df{9^  ±fi„)/d9=  0 exhibits 


a negative  wall  pressure,  P|,<0.  For  this  reason, 
the  onset  of  separation  was  determined  from  the 
condition  P( 9=  ±e„) 3 0.’  Similarly,  here  we  have 
associated  the  onset  of  separation  in  flows  with 
stress  relaxation  with  the  vanishing  of  the  wall 
pressure,  Pn'O.  The  observed  stabilizing  effect 
of  the  magnetic  field  is  due  to  the  increase  in  the 
wall  pressure  (Fig.  6)  with  Increasing  Hartmann 
number  X'.  In  conventional  incompressible  fluid 
dynamics  without  stress  relaxation,  the  conditions 
df{9=  ±9g)/d9’=  0 and  P(9*  0 lead  to  the  same 

separation  criterion.  j 

Interest  in  this  theoretical  problem  arose  in  j 

connection  with  experiments  on  boundary-layer  j 

separation  in  incompressible  liquid  metal  flow  | 

and  subsonic  magnetogasdynamic  flow  in  nonuni-  I 

form  magnetic  fields  and  ducts.’’’  If  the  Hart-  j 

mann  number  is  set  to  zero,  the  closed-form  solu-  | 

tions  presented  reduce  to  those  for  the  flow  of 
electrically  nonconducting  ordinary  fluids  with 
viscous-stress  relaxation  in  diverging  ducts. 
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APPENDIX  VTStOLS  STRESS  TRANSPORT  EOGATION 

Incompressible  magnetohydrodynamics  is  appli- 
cable to  conducting  liquids  such  as  liquid  metals, 
and  also,  as  an  approximation,  to  collision-domi- 
nated ionized  gases  and  plasmas  at  subsonic  flow 
speeds.  In  each  case,  the  viscous  momentum 
transport  is  due  to  the  heavy  atomic  particles 
while  the  electrons  affect  only  the  stress  relaxa- 
tion time  r through  cross-collisions.  For  this 
reason,  conducting  fluids  can  be  described  by 
one-fluid  magnetohydrodynamic  equations. 

The  stress -relaxation  equation  for  incompres- 


FIG.  5.  Critical  duct 
angle  Pj  for  separation  vs 
P for  various  3C,  with  ( — ) 
and  without  ( — ) viscous- 
stress  relaxation. 
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relaxation  equation  can  be  derived  from  the  Bom- 
Green  kinetic  equation  which  differs  from  the 
Boltzmann  equation  through  the  collision  integral 
for  many-body  collisions.  The  most  elegant  meth- 
od of  solution  for  the  Boltzmann  equation  ia  the 
13-moment  approximation  due  to  Grad/"*"  which 
Is  mathematically  also  more  rigorous  than  the 
Chapman-Enskog  and  Hilbert  methods."  The 
closure  of  the  13-moment  approximation  is  forced 
by  truncating  the  third-order  heat-flux  tensor 
In  terms  of  the  heat-flux  vector  ^ 


//£■ 


= f(<7,6rt+<?^6,i  + (7,8u£ 


(Al) 


FIG.  6.  Wall  pressure  Pj  vs  K for  various  Pieo'S”). 


slble  (V  0)  flows  is  derived  as  a special  case 
of  the  stress -relaxation  equation  for  compressible 
flows  based  on  the  Boltzmann  equation 
presuming  that  the  fluids  are  sufficiently  dense  and 
collision  dominated.  For  liquids  a similar  stress 


where /(c,r,f)  Is  the  distribution  of  thermal  vel- 
ocities c.  This  truncation  affects  mainly  the 
heat-flow  dynamics,  and  is,  therefore,  an  excel- 
lent approximation  for  quasi-isothermal  magneto- 
hydrodynamic flows  (q^o). 

Following  the  original  deduction  of  Grad  for  a 
neutral  one -component  gas,"  we  multiply  the 
Boltzmann  equation  for  particles  of  mass  m,  and 
charge  in  an  electromagnetic  field  E - B by 
imjc^c, and  integrate  it  over  the  entire 
space  / / /d*c,  of  the  thermal  velocities  c,  of  the 
r particles.  Thus,  the  following  moment-conser- 
vation equation  is  obtained  for  the  nonhydrostatic 
stress  tensor  n,«P,  -p^S  of  the  r component 
(summation  s over  the  remaining  components  s of 
the  fluid); 


^n,+  n,  X W,+  [n,  x + v,  • vn,+  n,v . . v,6)+  |{v5,+  [v5,]-‘  - fv  .5,6) 


+ (ff, . [n, . VvJ-  - !n, : v7,6)  jj - i5j6]C,.dS„  (A2) 


where 


C,.-  /•••  f if:f:-frf.K.igr.,^)gr.dndS.  (A3) 

is  the  binary  collision  integral,"  and  w,- 
Is  the  gyration  frequency  ([  ]**  designates  the  in- 
verse tensor).  In  the  13-moment  approximation, 
the  distribution  functions  /,(?„?,!)  and  /,(?„?,/) 
are  expanded  in  Hermlte  tensorial  polynomials, 
the  expansion  coefficients  being  the  first  (scalar) 
13  moments  of  the  distribution  function."  By 
means  of  these  expansions  and  Eq.  (A3)  it  can  be 
shown  that  Eq.  (A2)  is  of  the  form 


sn, 

SI 


(A4) 


where  is  the  viscous-stress  relaxation  time" 
and  T,.,  is  the  relaxation  time  describing  the  linear 
momentum  exchange"  between  the  components  r 
and  s*r.  6,.  is  an  abbreviation  for  the  remaining 
tensor  terms  on  the  left  side  of  Eq.  (A2),  and 
a,.,  are  numerical  coefficients.  For  representative 
times  t»r^,  the  term  STl^/st  in  Eq.  (A4)  is  negli- 
gible, and  one  obtains  the  qiuisi -equilibrium,  I),, 
proportional  to  the  sum  of  the  various  driving 
force  tensors."  The  tensor  §,  reduces  not  always 
to  the  velocity  gradients  ^v,  of  the  Navier-Stokes 
relation,"  and  for  representative  times  the 

13-moment  approximation  does  not  approximate 
the  phenomenological  Navier-Stokes  relation  and 
gives  better  results." 

For  magnetohydrodynamic  applications,  a sim- 


33 


...-V  . — . 


j pie-stress  transport  equation  (or  the  electrically  of  the  inverse  relaxation  times  The 

I conducting,  incompressible  fluid  as  a whole  can  Reynolds  number  in  tne  preceding  sir.iilarity  ana- 

f be  deduced  from  Eq.  (A2)  by  neglecting  those  lysis  and  t are  related  by  fl(0)«  [pi/(0)r //>(0)|r-‘. 

I terms  which  are  small  compared  with  the  leading  Equation  (a5)  can  also  be  derived  directly  from 

(terms.  The  contribution  of  the  electrons  (e)  to  the  elementary  physical  arguments.  The  Navier- 

fluld  of  ions  (i)  and  atoms  (a)  as  a whole  is  Insig-  Stokes  driving  force  -/>(Vv+  [VvJ"‘)  of  fl  follows 

nlficant  since  m,«  m,  , and‘*  1.  The  from  the  symmetry  argument  of  Einstein.  The 

magnetic  anisotropy  terms  are  negligible  since  term  aU  '3/  results  from  the  fact  that  the  Navier- 

u<,T,  v<  1 for  the  heavy  ions,  and  vanish  for  the  Stokes  quasi -equilibrium  n/T=  -p(rv  + [Vv]*‘)  de- 
neutral atoms  (u),  = 0).  The  V-v,  terms  can  be  velops  within  a time  of  the  order  of  the  “collision” 

disregarded  for  incompressible  fluids  and  sub-  time  t.  Finally,  the  convective  term  v*^fi  has 

, sonic  (compressible)  flows.  Theq,  terms  are  to  be  added  in  order  to  make  Eq.  (A5)  invariant 

j negligible  (or  quasi-isothermal  flows,  and  the  against  Galilei  transformations  (?'  = ?- wf,/'=/). 

; terms  n,  • Vv,  are  of  the  order  of  magrUtude  of  We  have  based  the  deduction  of  the  stress  trans- 
quadratic terms  in  Vv,,  and  therefore  small  com-  port  equation  on  Grad’s  13-moment  theory,  which 

pared  to  the  linear  ones.  The  stresses  in  the  elec-  gives  all  driving  forces  (for  the  vi.scous  stresses) 

tron  gas  have  no  effect  on  the  stress  distribution  which  have  a simple  physical  meaning,  i.e.,  not 

i of  file  fluid  as  a whole  (m,«m,_,).  Thus,  one  ob-  only  the  Navier-Stokes  forces  '^v  and  ~[Vv]*'. 

i tains  from  Eq.  (A2)  as  stress  transport  equation  It  is  seen  that  the  13-moment  theory  is  more  com- 

(or  incompressible,  quasi-isothermal  fluid  as  a prehensive  and  more  rigorous  than  the  Navier- 

i whole;  Stokes  theory,  and  “refuses  to  predict  results 

which  may  be  inaccurate.”  On  the  other  hand, 
9ll/9(+ V • vn  = -T"*n  -^(Vv+ [Vv]*‘),  (A5)  the  failures  of  the  Navier-Stokes  stress  equation 

may  assume  catastrophic  proportions,  e.g,,  it 

where  v«S«,m,v,/£n,w,  is  the  mean-mass  velo-  “predicts  smooth  solutions  for  shock  strengths  of 
city  of  the  fluid.  The  relaxation  frequency  t"'  of  infinite  magnitude  (with  a transition  from  negative 

the  total  stress  tensor  H is  a linear  combination  to  a positive  density).” 
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Self-similar  magnetohydrodynamic  diffuser  flows  with 
induced  magnetic  fields 

S.  H.  Choi  and  H.  E.  Wilhelm 

Department  of  Elettrtcal  Engineering,  Colorado  State  University,  Fort  Collins,  Colorado  80S2J 
(Received  22  February  1977;  RimI  mtnuKripl  received  II  July  1977) 

The  needy  dilTutcr  Row  of  i viscoue  incompretiible  Duid  with  a finite  electneal  conductivity  in  the 
preaence  of  an  external  magnetic  field  ii  analyzed  by  meant  of  ximilarity  theory.  Power  series  expansions 
for  the  velocity  and  induced  magnetic  fields  are  presented,  and  discussed  as  dependent  on  the  Reynolds, 
Hartmann,  magnetic  Reynolds  numbers,  and  duct  angle  Short-circuit  current  and  the  critical  values  of 
magnetic  Reynolds  number  and  duct  angle  for  flow  separation  are  calculated  It  is  shown  that  the  induced 
magnetic  field  changes  the  velocity  profile  and  coniiderably  reduces  the  critical  duct  angle  at  which  flow 
separation  occurs.  The  critical  magnetic  Reynolds  number  fur  flow  separation  becomes  smaller  as  the 
Reynolds  number  increases. 


i 

I 


I 


I.  INTRODUCTION 

Extensions  of  the  Jeffery-Hamel  flow  theory'’*  to 
magnetohydrodynamics  have  been  made  by  a few  au- 
thors.'"* These  approaches  to  magnetohydrodynamic 
flows  in  a diffuser  are  based  on  the  approximation,  (i) 
that  the  induced  magnetic  field  is  negligible  (magnetic 
Reynolds  number  Af « 1),  ’"*  or  (li)  that  the  working  fluid 
la  an  ideal  electric  conductor  (electrical  conductivity 
o- 

In  the  following,  the  steady  flow  of  an  incompressible 
viscous  fluid  with  finite  conductivity  between  two  inclined 
walls  (diffuser)  in  the  presence  of  azimuthal,  external, 
and  Induced  magnetic  fields  is  analyzed.  The  electrical 
conductivity  is  assumed  to  be  constant  throughout  the 
diffuser.  This  approximation  is  applicable  to  conducting 
liquid  and  incompressible  (subsonic)  gaseous  plasma 
flows,’ 

It  is  shown  that  the  steady  motion  of  the  conducting 
fluid  does  not  induce  a magnetic  field  in  the  azimuthal 
direction,  but  does  induce  it  in  the  radial  direction.  An 
analytical  solution  to  the  magnetohydrodynamic  flow 
problem  is  accomplished  by  reducing  the  nonlinear  par- 
tial differential  equations  to  ordinary  ones  via  a similari- 
ty transformation,  and  solving  the  latter  equations  in 
terms  of  power  series.  The  magnetohydrodynamic 
field  profiles  are  presented  in  dependence  of  the  Rey- 
nolds (R),  Hartmann  (H),  magnetic  Reynolds  (M)  num- 
bers, and  duct  angle  (^g). 

In  the  limiting  case  of  small  magnetic  Reynolds  num- 
bers (Af « 1),  the  nonlinear  boundary  value  problem  re- 
duces to  the  conventional  magnetohydrodynamic  diffuser 
problem,*  which  is  based  on  the  approximation  that  the 
induced  magnetic  field  is  small  compared  with  the  ap- 
plied magnetic  field. 

II.  PHYSICAL  PRINCIPLES 

The  electrically  conducting  fluid  flows  between  two 
nonparallel,  insulating  walls  and  two  parallel  electrode 
plates  as  shown  in  Fig.  1.  The  electrodes  are  short- 
circuited  so  that  there  is  no  potential  difference  between 
them.  The  influence  of  a load  in  the  external  circuit  has 


been  investigated  previously.*  The  conducting  fluid  is 
injected  through  the  inner  duct  section  at  r = Vg  and  re- 
moved downstream  through  the  outer  duct  section  at  r 
= ri.  The  boundary  layers  at  the  electrodes  are  disre- 
garded compared  with  those  at  the  insulating  walls  pre- 
suming that  the  inter-electrode  spacing  is  large,  r. 

» 2(rg  + >'i)eo  . The  external  magnetic  field  has  its  source 
in  an  electric  current  /g  flowing  through  an  infinitely  long 
conducting  rod  located  in  the  z axis.  In  accordance  with 
Ampere’s  law,  • ffs  = /.ig /g , the  external  magnetid  field 
is  azimuthal  and  has  the  induction  (pg  is  the  permeabili- 
ty of  vacuum), 

iig  = (fig/2ir)(/g/>-)  . I 

In  the  absence  of  flow  sources  or  sinks  at  the  walls, 
the  velocity  field  is  radial,  i.e.,  v=  [i;(»-,  6),  0,  O].  The 
flow  of  the  fluid  across  the  external  magnetic  field  Bg 
induces  a current  density  field  in  the  axial  direction,  1 

J,  = (jiE,  + vBg)  presuming  that  the  Hall  effect  is  negligible.  | 

A magnetic  field  is  induced  in  the  radial  direction  by  j 

the  current  density  if, . Since  V “E  = 0 aixi  V • J = . (E  | 

+ v»B)  = a^-E  = 0,  a homogeneous  electric  field  E is  pos-  j 

sible.  However,  E,  = 0,  since  the  terminal  voltage  is  | 

zero  due  to  the  short-circuited  external  circuit,  and  | 

=0  since  the  boundary  conditions  for  £,,,  are  ho-  i 

mogeneous.  The  Hall  effect  is  assumed  to  be  negligible,  j 

uiT«  1 {w=eB/in  is  the  gyration  frequency  and  t is  the  1 

collision  time  of  the  electrons).  Thus,  Maxwell’s  equa-  ] 
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tions  and  Ohm's  law  give,  for  the  magnetic  field  B = , 

fif  , 0),  the  induction  equation, 

7*8  . 


III.  BOUNDARY  VALUE  PROBLEM 


In  nundimensional  form,  indicated  by  the  overbars, 
the  radial  velocity  F(f,  6),  the  pressure  p{r,  9),  and  the 
magnetic  field  fi),  B,(r,  »),  0]  are  described 

by  the  nonlinear  boundary- value  problem; 


1 ~ (?r)»0  , 

r 9r 


_ 9v  86  1 r 9*r  1 9P  1 


8f 


8r  R[.9r*  r 9r  r* 


d*t7 


H*  1 fl  9 - . 1 8Jf,1- 


]b,. 


(1) 


(2) 


i 

? ~»e  * 


2 J_  8r  «•  1 1 r » 1 9B,1- 

R F*  Be*  R Mr  L8F  ~ 86  J ’ 


J-fi  -L 

Af  9F 


1 9 /--ri  1 1 

■=  T=  {rBr)+- 
r Br  r 


■]  =l’B*  . 


BB, 

Be 


= 0 . 


(3) 

(4) 

(5) 


In  Eqs.  (l)-(5),  the  Reynolds  (A),  Hartmann  (A),  and 
magnetic  Reynolds  (Af)  numbers,  and  the  dimensionless 
variables  are  defined  by 

"*  = (?&  ' Af  = awa.  (6) 

and 

F = i'A’o,  P=p/!Hv\,  B,  = B,/(tigIo/2irro)  , 

= /(po/(»/2irro)  , r = rfr^,  (7) 

where  vq  is  the  input  velocity  at  r=ra,  and  6=0  (cen- 
tral stream  line).  The  solutions  to  Eqs.  (l)-(5)  are 
subject  to  the  following  conditions: 


v(r,  ± 6o)  = 0 , 

(8) 

B,(r,  ± Bg)  = Bg(r^  , 

(9) 

v(r,  +e)  = v(r,  -e)so  , 

(10) 

A,(F,  +6)  = -A,(F,  -6)  , 

(11) 

S-.f'^^-rde. 

(12) 

PgVgYg 

Equation  (8)  contains  the  boundary  conditions  for  no  slip 
of  v{r,  e)  at  the  walls.  Equation  (9)  considers  the  con- 
tinuity of  the  normal  component  B,  at  the  wall  inter- 
faces (A,(F)  = 1/F).  Equations  (10)  and  (11)  represent 
the  symmetry_condition  lor  i7(F,  6)  and  the  asymmetry 
condition  for  B,{r,  6),  respectively  [Eq.  (11)  will  be  jus- 
tified later].  Equation  (12)  gives  the  flow  rate  Q per  unit 
depth  (Ae  = 1). 


IV.  SIMILARITY  TRANSFORMATION  AND  SERIES 
EXPANSIONS 

The  continuity  equation  (1)  indicates  that  the  velocity 
field  is  of  the  self- similar  form 


1822  Phv$.  Fluids,  Vol.  20,  No.  11,  November  1977 


36 


r(F,  fl)=f(e)/F . (13) 

For  dimensional  reasons,  the  pressure  and  magnetic  in- 
duction fields  are  written  as 

p(F,  e)=G(e)/F‘ , (14) 

B,(r,  6)  = ip(e)/r  , (15) 

B,(F,  fl)  = (l/F)[l + 0(6)1  . (16) 

In  Eq.  (16),  0/F  represents  the  induced  magnetic  field 
in  the  aslmuthal  direction,  and  1/F  is  the  external  field. 

Substitution  of  Eqs.  (IS)  and  (16)  into  Eq.  (5)  (con- 
tinuity equation  for  B)  gives 

d(t>/de  =0  . (17) 


Hence, 

0(6)  = const5O,  andB,  = l/F  (18) 

by  Eq.  (9).  Thus,  the  azimuthal  component  of  magnetic 
induction  is  not  affected  by  the  steady  motion  of  the  con- 
ducting fluid.  Substitution  of  Eqs.  (13),  (14),  (15),  and 
(18)  into  Eqs.  (l)-(5)  yields  the  coupled  nonlinear  ordi- 
nary differential  equations: 


(19) 

(20) 

M de  ’ 

(21) 

where 

/•(6=±6o)=0,  F(+e)=F(-fi)so  , 

(22) 

0(+6)  = -0(-6)  , 

(23) 

and 


Fde=  ^ . (24) 

R n 

Since  the  velocity  field  is  symmetrical  and  greater  than 
zero  for  pure  outflow  [/■(+6)  = F(-  6)^0],  di),/dp  s 0 
everywhere  and  0(+ 6)  = - 0(- 6)  by  Eq.  (21).  Thus,  the 
condition  in  Eq.  (23),  or  Eq.  (11)  is  justified. 

Integration  of  Eq.  (20)  with  respect  to  6 yields 

where  gg  is  an  integration  constant.  In  the  limit  of  large 
Reynolds  number  A,  G(6)  in  Eq.  (25)  becomes  the  con- 
stant i.e.,  the  pressure  field  would  no  longer  depend 
on  the  transverse  coordinate  (6)  (in  agreement  with  the 
so-called  boundary-layer  approximation'). 

Elimination  of  G(6)  in  Eq.  (19)  by  means  of  Eq.  (25) 
leads  to 


+AF*+4F+ 


M 


+ Co  = 0 


(26) 


where  Cg  ^ 2Rgg  is  an  integration  constant  proportional 
to  go  • Hence,  the  problem  reduces  to  solving  two  cou- 
pled nonlinear  differential  equations  [Eqs.  (21)  and  (26)]. 


In  the  limiting  case  of  vanishing  induced  magnetic 
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TABLE  1,  IntCKratlon  constant  Co  = 2i{  go 
forR-10*,  .M  = l,  60=5°,  and  various  H. 


Co 

•IR 

l.OlOl  X 10< 

3R 

2.0004  X 10* 

6R 

3.  9999  X 10* 

lOR 

8. 9999  X 10* 

Substitution  ol  Eq.  (35)  into  Eq.  (28)  pves 
as=i[Co  + flai-(4-H*)a,l  , 
nj=^[6/eoirtj-3(4-ff‘)a,-Afff*afl  . < 

+9a*)-  5{4- «*)«,- 2AlH‘n,<i,l  , 
««  + 3CVijn,)  - 7(4  -H*)a, 

-Aftf*(2a,a5+oJ)l  , 


Meld  (M  -0),  by  eliminating  dtli/de  in  Eq.  (26)  by  means 
of  Eq,  (21)  and  taking  i(i(e)-0  and /■(e)-/(0)//i,  Eq.  (26) 
reduces  to  the  equation,  /"  + (4 *)/+/* +a  =0,  which 
has  been  solved  in  closed  form  In  terms  of  elliptic  func* 
tlons.^ 

Normalizing  1/1(6)  by  the  factor  Af, 

<li(e)=-ip(0)/M  , (27) 

and  eliminating  /"(fl)  from  Eq.  (26)  by  means  of  Eq.  (21), 
one  obtains  a single  nonlinear  inhomogeneous  differential 
equation  for  4'(6) 

^ = . (28) 

By  Eqs.  (21)  and  (27),  the  conditi  ns  in  Eq.  (22)  and  (23) 
reduce  to 

'I-(0)  = 0,  (29) 

dM±6o)/de  = 0 , (30) 

£f*4'(0)/rfe*  = 0 . (31) 

The  flow  rate  is  given  by 

0 r'o  _ f •♦o  /d*\  ..  ..... 


-s.  = r 


C **0  /d*  \ 

J..0  \de) 


de=T2'K±8o) 


l.e., 

Q = ^2jiRM±ea)  . (32) 

Equation  (32)  could  be  used  as  a boundary  condition  for 
'l'(6)  presuming  that  the  flow  rate  Q is  given  (R  would 
then  no  longer  be  independent).  It  is,  however,  mathe- 
matically more  convenient  to  assume  the  Reynolds  num- 
ber R(0)  of  the  central  streamline  to  be  given  by’ 

R^ R{0)=p^v{r,  B =0)r/^l  , (33) 

which  implies 

F{B  =0)  = 1.  (34) 

The  solution  <i'(6)  of  Eq.  (28)  can  be  expanded  in  terms 
of  a power  series  with  respect  to  6,  which  satisfies  the 
conditions  in  Eqs.  (29)  and  (31) 

♦(e)  = Z]  ai«i  e*"*‘ . (35) 

mO 

The  equivalent  condition  to  Eq.  (34)  is,  by  Eqs.  (21)  and 
(27), 

d*/de  = - 1 at  e = 0 . (36) 

Equation  (36)  determines  the  first  power  coefficient  in 
Eq.  (35), 

a,=-l,  (37) 


° ^J«-T)  - 1 )(4  - 

^■R  +l){2n- 2k- 

*•0 

-1 

1 • '•■1>2,  3, 

*■’  ■’  (39) 

where  S,„  = 1 for  n = 1 and  fli,  = 0 for  n # 1 . Equation  (39) 
is  a two-term  recurrence  relation  from  which  (or  a 
given  one  can  readily  compute  o,„^  and  then  a,„j , 

^ Ol'.  desired. 

The  nondimensional  functions  ♦(e)  and  f(B)  =-  ♦'(e) 
are,  by  Eqs.  (35)  and  (37), 

♦(e)=-e+ ^ a,,„,e*"*‘ , (40) 

nal 

/'(e)  = l-  E (2«  + i)a,„..e‘" . (41) 

n»i 

The  condition  in  Eq.  (30)  determines  the  integration  con- 
stant Cg  in  Eq.  (39), 

1=^  (2n  + l)a2„,,  ej"  , nt,»i  =ni,Hi(^o^  • (42) 

fial 

Equation  (42)  is  a function  of  C,  for  a given  duct  angle 
6q  . By  Eqs.  (25)  and  (27),  the  pressure  field  G(e)  is, 
in  terms  of  ♦(e)  given  by  (Cci*2Rffo). 

G(e)  = [2f(e)-Afi/’^’(e)/2  + Co/2)  . (43) 

Equations  (39)- (43)  represent  the  formal  mathematical 
resolution  of  the  boundary- value  problem  under  con- 
sideration. 


V.  DISCUSSION  AND  RESULTS 

Based  on  a numerical  iteration  of  Eq.  (42),  the  values 
of  the  integration  constant  Cg  are  given  for  various  Rey- 
nolds (R),  Hartmann  (R),  and  magnetic  Reynolds  (Af) 
numbers  in  Tables  I and  11. 


A.  Flow  fisWt 

In  Fig.  2,  the  nondimensional  velocity  amplitude  F(B) 
is  shown  for  R =10*,  .Vf  = l,  eg  = 5°,  and  various  values  of 
Af.  The  width  of  the  boundary  layer  is  considerably  re- 
duced as  the  Hartmann  number  R Increases.  The  in- 
fluence of  the  magnetic  Reynolds  number  .Vf  on  the  ve- 
locity profiles  is  exhibited  for  R •10*,  10’,  and  R*  = 3R, 
eg  = 5°  in  Figs.  3 and  4,  respectively.  It  is  seen  that 
the  velocity  distribution  becomes  flatter  in  the  center, 
but  the  boundary-layer  width  remains  almost  constant 
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TABLE  II.  Integration  conatant  C^rn'iHgf^  for  =3R,  0^  ^ 5°, 
and  various  M 


R 

M 

to* 

10* 

10* 

1 

2.1226  X 10* 

2.0004  X 10* 

2.0000  X 10* 

to 

2.1416  X 10* 

2. 0055  X 10* 

2.0006  X 10* 

2U 

2 1620  X 10* 

2.0107  X 10* 

2.0012  X 10* 

24.6 

30 

2.  1810  X 10* 

2.0167  X 10* 

2.0014  X 10* 

42.6 

50 

222. H7 

2.2201  X 10* 

2.4087  X 10* 

2.0215  X 10* 

as  A/  decreases.  The  velocity  profiles  for  the  critical 
magnetic  Reynolds  numbers,  above  which  flow  separation 
occurs,  are  plotted  with  a dotted  line  in  Figs.  3 and  4. 

In  Figs.  5 and  6,  the  nondlmensional  induced  magnetic 
field  amplitude  4’(6)[>  (/>(e)/A/]  is  shown  for  R = 10’,  10^ 
with  H*  = ZR  and  60  = 5°,  respectively.  In  these  figures, 
4>(0)  increases  in  magnitude  from  the  center  to  the  walls 
for  Af  less  than  a certain  critical  value,  for  which 
flow  separation  sets  in.  'i(e)  for  the  critical  magnetic 
Reynolds  number  is  plotted  with  dotted  line.  Unless 
A/  is  of  the  order  of  unity  or  larger,  it  is  seen  that 

is  small  compared  with  the  applied  field.  In  pure 
outflows,  the  induced  magnetic  field  does  not  vanish  at 
the  walls  since  the  net  current  flow  across  the  diffuser 
in  the  axial  direction  is  nonzero. 

In  accordance  with  Eqs.  (6),  (7),  (13),  and  (18),  the 
axial  current  density  is 

= (44) 


FIG.  3.  F(9)  versuB  eforR  = 10*,  lP  = 3R,  #0  = 5",  and  vari- 
ous M.  J,  is  proportional  toF(0l. 


9/80 


FIG.  4.  F(0)  versus  e for/?  = 10*.  lf'  = 3R  0(1  = 5”,  and  vari- 
ous M.  J,  is  proportional  to  F(9). 


in  dimensional  form,  or  in  dimensionless  form 


(45) 


H. 


ous  M.  la  equal  to  Af  I 4'  (4  Oo)  I . 


1824 
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FIQ.  6.  _*(«)  versus  «for  = = »o  = 5°.  and  vari- 

ous M.  /„  la  equal  to  4f  I 4'  (i  fg)  I . 

or  in  terms  of  the  flow  rate, 

7„  = {M/R)Q  ln(ri  / rg)  , (47) 

where  7„  =/„  /(/„ /27r)  and  Q = Q/n  =R  f:Sojr(0)de.  In  ac- 
cordance with  Eq.  (32),  the  flow  rate  is  given  in  terms 
Of  *(*60). 

/,c  = M I ♦(*  eg)  I ln(n  /rg)*  . (48) 

Thus,  the  axial  current  density  in  Eq.  (45)  and  net  short- 
circuit  current  in  Eq.  (48)  can  be  visualized  for  various 
R,  H,  and  M numbers  with  eg  = 5°  in  Figs.  3,  4 and  5,  6, 
respectively.  An  additional  illustration  for  is  given 
in  Fig.  7 as  a function  of  M for  various  values  of  H. 

B.  Flow  isparition 

Magnetohydrodynamic  flow  separation  sets  in  if  one  of 
the  flow  parameters  goes  beyond  its  "critical"  value, 
e,  g. , if  the  duct  angle  fig  is  larger  than  a critical  value, 
0c  = 6g{R,  H,  M),  If  R,  H,  and  M are  assumed  to  be  giv- 
en, 6c  cannot  be  chosen  independently  since  it  assumes 
the  role  of  an  eigenvalue. 


FIG.  7.  ln(ri/rg)*|  versus  M for  various  H. 
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FIG.  8.  The  critical  duct  angle  flt>vr  separation  versus 

R for  H = 3 and  various  M,  with  ( ) and  without  ( ) 

Induced  magnetic  Induction. 


The  power  series  solution  derived  in  Eqs.  (40),  (41). 
and  (43)  describes  physical  flows  as  long  as  =C(6 
= ±eg)^0  for  any  set  of  parameters  R,  H,  M,  and  6g  , 
since  the  pressure  field  has  to  be  positive  everywhere, 
i.e.,  G(e)2:0,  lei^sfig.  Since  F(e)=0  at  the  walls,  G„ 
in  Eq.  (43),  which  is  related  to  the  wall  pressure,  is 
given  by 

G.  = (Cg-A/^**|.^^)/2fi  . (49) 

As  the  pressure  field  drops  to  zero  at  the  walls,  the 
laminar  flow  can  no  longer  be  realized,  so  that  flow 
separation  occurs  for  duct  angles  larger  than  a critical 
duct  angle  9^  at  which  separation  first  sets  in.*  Hence, 
at  the  onset  of  separation,  Eq.  (49)  yields 

Cg  = ^ . (50) 

For  the  limiting  flow  of  vanishing  wall  pressure,  Eq. 
(50)  reduces  Eqs.  (40)  and  (41)  to 

-(Cg/Afi/*)‘'‘=-e.+  X;  , (51) 

ff*l 

and 

1 = ^ (2n +l)o,„.,eJ"  , (52) 

where  =a,,.i(Cg).  The  critical  duct  angle  9^  =9,(71, 
M,  H)  is,  therefore,  obtained  by  solving  the  two  coupled 
equations,  (51)  and  (52),  numerically.  In  Fig.  8,  the 
critical  duct  angles  are  plotted  as  a function  of  R for 
various  values  of  Af  and  a given  H.  It  is  seen  that  9^ 
decreases  as  R increases.  As  M increases,  the  critical 
duct  angle  9,  decreases  for  low  Reynolds  numbers  {R 
5 10)  and  increases  for  intermediate  Reynolds  numbers 
(I0<il  <10’),  while  9j  varies  little  with  Af  > 1 at  high 
Reynolds  numbers  (H  >10*).  The  critical  duct  angle  9, 
for  the  limiting  case  of  vanishing  induced  magnetic  field 
(Af-  0)*  is  shown  with  a dashed  line  in  Fig.  8.  Compari- 
son indicates  that  the  induced  magnetic  field  reduces 
the  critical  duct  angle  by  about  50%. 

A simple  way  of  simulating  the  self- similar  flow  in- 
jection is  to  put  the  diverging  channel  with  free  entrance 
and  exit  openings  into  a larger,  similar  duct  through 

S.  H.  Choi  and  H.  E.  Wilhelm  1825 


which  the  conducting  fluid  is  pumped.  This  method  is 
commonly  used  in  experiments  on  ordinary,  noncon- 
ducting Kamel  flows  between  inclined  walls. 
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Stark  effect  in  three-dimensional  stochastic  electric  fields 
and  a static  magnetic  field  of  a helium  plasma*’ 
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Dt^rtment  ^ Suclmr  VtUnnU^,  Sum/  Hi.  Kona 

(Rscslvcd  14  December  1974;  accepted  for  puUicslion  II  Apnl  1977) 

Atomic  abaorptKNi  epectra  of  the  He  I plaama  in  the  Matic  magnetic  Held  in  the  absence  of  laser  radiation 
are  different  from  spectra  in  the  presence  of  laser  radiation.  Especially,  it  has  been  theoretically  predicted 
that  the  spectral  shift  of  the  satellite  line,  which  is  seen  in  the  latter  case,  does  not  exist  in  the  former 
CMC  In  this  paper,  the  effect  of  ihree-dimcnsMNkst  hi^-frequency  electric  fields  on  the  cr^,  a. . and  ir 
components  of  the  allowed  1 'So-3  'P|  transition  and  the  forbidden  1 ’5o‘3  ‘Dj  lines  in  atomic  absorption 
spectra  of  the  He  I plasma  m a static  magnetic  field  in  the  absence  of  a laser  beam  is  investigated  from 
the  coupled  equations  for  the  probability  amplitudes  of  atomic  states.  It  is  shown  that  they  should  give 
riac  to  a spectral  shift  of  the  allowed  lines  and  satellites  about  forbidden  lines.  The  calculaticm  reveals  that 
the  shift  is  proportional  to  the  mcan*squarc  stochastic  electric  field,  snd  that  the  ratio  of  the  intensity  of 
the  near  satellite  to  that  of  the  far  satellite  is  different  from  the  ratio  obtained  by  the  second-order 
perturbatioo  theory  or  the  three-level  rotating-wave  approximation.  Discussion  is  given  for  the  validity  of 
this  theory. 

PACS  numbers:  3160. -t-i.  32.70Jz.  92.70.Kz 
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I.  INTRODUCTION 

The  tranaitlons  betweea  the  bound  states  of  an  atomic 
species  in  the  coUlstonless  plasma  are  mainly  produced 
by  absorbing  or  emitting  an  optical  photon  of  the  radia- 
tion field  or  the  multiple  microwave  photons  of  the 
turbulent  electric  fields. 

If  the  radiation  field  is  composed  of  only  the  back- 
ground radiation  field  which  covers  a wide  spread  of 
frequencies  and  polarizations,  all  upper  levels  of  the 
allowed  transitions  are  resonantly  pumped  by  the  back- 
ground radiation  field.  Accordingly,  the  probability  for 
finding  the  atom  in  the  upper  level  of  the  allowed  transi- 
tion is  greater  than  those  In  the  upper  levels  of  the  for- 
bidden transitions.  Hence,  the  probability  for  finding 
the  atom  in  the  upper  level  of  the  forbidden  transition  is 
strongly  dependent  on  those  In  the  upper  levels  of  the 
allowed  transitions  which  are  coupled  with  the  former 
level  by  turbulent  electric  fields  according  to  the  elec- 
tric dipole  selection  rules.  In  this  case,  the  turbulent 
electric  fields  produce  primarily  Stark  shifts  of  the 
allowed  lines  but  do  not  significantly  shift  those  of  the 
forbidden  lines  ("asymmetrical  Stark  repulsion"). 

Theoretical  treatments  of  the  high-frequency  Stark 
effect,  based  on  second-order  time- dependent  perturba- 
tion theory,  have  been  given  by  Baranger  and  Mozer, ' 
Relnhelmer,  * and  Copper  and  Rlngler.  ’ Kunze  et  al.  * 
eTctended  the  perturbation  calculation  up  to  fourth  order. 
However,  these  theories  did  not  Include  Stark  shifts  of 
the  levels.  Recognizing  this,  Hicks  ef  al. ' have  devel- 
oped a theory  which  can  include  Stark  shifts  (actually 
profile).  Their  theory  is  based  on  two  consecutive 
perturbation  calculations,  the  first  calculation  being  for 
the  electric  field  and  the  second  calculation  being  for 
the  zero  radiation  field  (spontaneous  emission).  They 
have  neglected  the  fact  that  at  least  one  intermediate 
upper  state  Is  not  only  coupled  with  the  initial  upper 
state  by  the  interaction  with  the  electric  field  but  also 
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with  the  final  lower  state  by  the  induced  absorption  or 
the  induced  and  spontaneous  emission  of  a photon. 
Therefore,  their  theory  is,  in  principle,  valid  only  in 
the  absence  of  the  strong  radiation  field.  Furthermore, 
they  assumed  that  the  upper  states  have  equal  popula- 
tions, a fact  that  is  certainly  not  in  general  true. 

Recently,  Prosnitz  et  al.‘  proposed  a (nonperturba- 
tioh)  theory  which  eiqilains  rather  well  the  line  shape  of 
both  the  satellite  and  allowed  line  in  the  presence  of 
laser  radiation.  Their  theory  is  based  on  the  coupled 
equations  for  the  probability  amplitudes  with  phenome- 
nological damping  terms,  but  they  do  not  Include  any 
term  for  the  optical  excitation  or  the  radiation  field.  In 
their  theory,  the  spectrum  of  the  satellite  lines  would 
be  the  same  whether  the  laser  radiation  is  present  or 
not.  Hence,  it  is  obvious  that  their  theory  cannot  explain 
most  characteristic  effects  caused  by  the  radiation 
field.  In  particular,  their  theory  will  give  rise  to  in- 
valid results  for  the  spectral  positions  and  intensities  of 
the  satellite  lines  in  the  absence  of  the  laser  field.  Even 
in  the  presence  of  the  laser  radiation,  the  satisfactory 
agreement  between  experiment*  and  theory*  appears,  in 
part,  to  be  due  to  the  use  of  adjustable  phenomenologi- 
cal parameters. 

For  these  reasons,  I recently  have  developed  a 
comprehensive  theory  which  is  valid  for  strong  electric 
fields  and  includes  Stark  shifts  of  the  levels  for  a three- 
level  atom  with  two  closely  spaced  upper  levels,  both 
in  the  absence  and  in  the  presence  of  the  laser 
radiation.  ’ 

The  formula  for  a spectral  shift  of  the  allowed  line  In 
this  theory  is  not  obtained  in  the  theories  mentioned 
above.  Further,  this  theory  explains  the  the  experimen- 
tal observation  that  the  ratio  of  the  intensity  of  the  near 
satellite  to  that  of  the  far  satellite  in  the  absence  of  a 
laser  beam  is  different  from  the  one  predicted  by  the 
second- order  perturbation  theory, ' or  the  three- level 
rotating-wave  approximation, ' and  that  the  Stark  repul- 
sion Is  asymmetrical. 
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This  paper  Is  an  extension  o(  the  previous  theory'  to 
Include  more  than  two  upper  levels  and  the  Interaction 
of  a mafnetlc  field  with  the  excited  atom.  I applied 
the  theory  to  the  o.,  o.,  and  v components  of  the  allowed 
1 '.‘<•-3  'P|  line  and  the  forbidden  1 '.S,-3  'D,  lines  of  the 
magnetic  He  I plasma  In  the  absence  of  laser  radiation. 
For  mathematical  convenience,  I treat  the  absorp- 
tion spectrum  rather  than  the  emission  spectrum.  The 
emission  spectrum  can  be  obtained  similarly  by  taking 
the  upper  state  as  the  initial  state.  Generally,  both 
spectral  profiles  are  not  the  same  since  the  coUlslon- 
less  turbulent  plasma  Is  not  in  local  thermodynamic 
equilibrium. ' 

II.  THEORY 

The  time-dependent  SchrOdlnger  equation  for  the 
helium  atom  and  the  radiation  field  with  the  Coulomb 
potential  V(r)  and  the  scalar  and  vector  potentials  e 
and  A Is 

+ ,*  + vr(^)  + ^ (V  * A) . 8 + 

xk(/)>.  (1) 

where  H,  la  the  quantized  tree  radiation  field  Hamil- 
tonian and  S la  the  spin  angular  momentum  operator. 

Here,  the  single-particle  (one-electron)  state  Is  used 
(or  describing  the  atomic  state  since,  in  moat  optlcal- 
(leld-produclng  and  turbulent-fleld-producing  transi- 
tions, the  initial  and  (Inal  states  differ  only  In  the  mo- 
tion at  one  electron.  Further,  the  central  field  approxl- 
matkii"  Is  used  (or  the  electrostatic  interactions  be- 
tween two  electrons  and  the  spln-orblt  Interaction  Is 
naglected  (Russall- Saunders  state). 

In  Eq.  (1),  the  scalar  and  vector  potential  can  be 
decorapoeed  In  the  form 

A = A*  + A*  + A*  (A*  = 0),  (2) 

(**  = **  = 0),  (3) 

where  superscripts  R,  M,  and  S denote  the  radiation, 
static  magnetic,  and  stochastic  electric  fields, 
respectively. 

Since  the  stochastic  field  E*  is  homogeneous  within 
atomic  dimensions.  Its  potential  * can  be  written 

♦(f)  = -r.B(/),  (4) 

where 

£y(f)=  S Fjcoe(fM  + ^^)  (/=x,y,*)  (5) 

are  the  vector  components  of  the  stochastic  field,  £j 
are  Its  amplitudes,  O are  its  frequencies,  and  are 
Us  random  phases.  For  convenience,  we  omit  the  sum- 
mation notation  In  Eq.  (5). 

The  vector  potential  A*  (or  a uniform  magnetic  field 


at  induction  B is 

A'  = iBxr, 

(6) 

V-A*  = 0, 

(7) 

(ieM/m  lA* . W = - («/2m  )B . L, 

(8) 

where  L=-  fSrxv  is  the  angular  momentum  operator. 
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frequencies  and  polarizations  with  no  phase  relations 
between  the  different  frequency  components  and  polari- 
zations. Since  it  is  Inconvenient  to  use  the  summation 
notations  (or  different  frequencies  and  polarizations 
everytime  we  write  the  radiation  field,  we  will  choose  a 
typical  monochromatic  and  linearly  polarized  field  to 
calculate  the  probability  of  a transition  between  station- 
ary states  and  omit  the  summation  notation  Itherefore, 


we  should  bear  in  mind  that  A(r)  is  the  summation  of 
such  fields].  The  vector  potential  A In  the  Schrodlnger 
picture  (or  the  monochromatic  and  linearly  polarized 
wave  field  is" 


A(r)  = (»/2f 0 TO))'  '^e[a  exp(ik . r)  + a*  exp(-  fk . r)],  (9) 


where  ci  is  the  dielectric  constant,  a and  o*  are  annihila- 
tion and  creation  operators,  respectively,  r Is  the  sys- 
tem volume,  (1)  is  the  frequency  of  the  wave  mode,  I Is 
the  photon  polarization  vector,  and  k is  the  propagation 
vector. 

The  term  H,  in  Eq.  (1)  Is  the  energy  of  the  quantized 
source-free  radiation  field  in  the  absence  of  the  atom, 

H,  = *u.(a*a  + i).  (10) 

The  energy  eigenstates  of  the  radiation  oscillator  of 
polarization  direction  e,  momentum  Sk,  and  energy  Ru 
satisfy  the  eigenvalue  equation 

*(pa*o|w>  = guiY|W>,  (11) 

where  N Is  the  occupation  number. 

The  numerical  value  of  the  quadratic  term  (eA’/lm) 
is  negligibly  small  even  in  more  than  the  second-order 
perturbation  calculation.  Hence,  we  neglect  it  in  the 
following.  Then,  by  choosing  the  Z direction  the  same 
as  the  B direction,  the  Hamiltonian  in  Eq.  (1)  Is  written 

H=H,  + H',  (12) 

«,  = - (*V2m)V*  + V(r) - (e/2m)B,-  (L,  + 2S.)  + 

(13) 

tf'  = - (c/m)()l/2{,Tw)‘  '*[a  exp(ik-  r)  + o*exp(-  fk-  r)] 
x(e.P)-cr  •£(().  (14) 

An  uilperturbed  eigenstate  of  Hf  may  be  written 

|a,.V)  = |o>  |Af>  (15) 

with  eigenvalue 

where  a refers  to  the  atom. 

Let  the  wave  function  I p(t)>  of  the  perturbed  wave 
equation  (Eq.  (1)]  be  expanded  In  terms  of  this  com- 
plete orthonormal  set  {ls,Af>}  by 

lm)=z:  C,^„(f)exp(-«->E,^„f)|a.Af),  (17) 

l•,X) 

By  means  of  the  orthogonality  property  of  lo,  Af>,  one 
obtains  the  differential  equation  tor  the  expansion  coef- 
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V 


( 


I 


I 

I 

i 

j 


flclents  (probability  amplitudes), 

X esp(t(ii  1^  |,<,  <r>0,  (1 8) 

where 

“ Ifc  «),  !••,  *■>  = •■'(£  Ifc  »>  - £ !,•,  »•>) 

= + (N-Ar')u..  (19) 


If  we  use  the  definition  of  creation  and  annihilation 
operators  and  the  orthogonality  relations,  we  obtain  by 
taking  only  the  electric  dipole  term  In  expansion  of 
expdk.r) 

<a,  N\H'\a',N')  = ~ Uui i..,(«/2(,tu)’ ^’(« . (a | r In') ) 
x((Al'  + l)'/»#,,^.,  + (Af')‘'’8»,*..,l 
-eE-<a|r|a')«»,,r,  (20) 

where 


<*>|.).I.*>  = 8''(£|.,-£|.-,).  (21) 

since  {a,N\H'\a',N')  does  not  involve  the  spins  and 
the  spin  functions  for  different  S’s  are  orthogonal,  the 
spin  state  does  not  change  in  the  transition.  According- 
ly, It  will  not  be  specified  in  the  following. 

By  /f,  In  Eq,  (13),  the  unperturbed  singlet  atom  state 
in  the  static  magnetic  field  Is  specified  by  quantum  num- 
bers n,  J,  and  M (L  = J for  the  singlet  state),  and  Its 
energy  Is 

£|n/,i»)  = £?n/)  + W*a(t,  (22) 

where  <i>j,  Is  the  Larmor  frequency  and  Is  the  en- 
ergy of  the  level  \h,J)  when  no  magnetic  field  Is  applied 
to  the  plasma.  Here  the  intensity  of  the  static  magnetic 
field  B Is  limited  to  the  range 


B<B 


• * 


(23) 


where  S,  is  a critical  value  above  which  the  I'aschen- 
Back  effect  becomes  Important. 

The  selection  rules  of  the  singlet  helium  states  are, 
because  of  and  Eq.  (20), 

M'  = M,  M±l,  (24) 

J'  = J±1  (25) 

for  the  transitions  both  by  the  radiation  field  and  by 
three-dimensional  stochastic  fields. 

In  this  paper,  the  transitions  from  the  n=;l  level  to 
the  >1  = 3 level  of  the  singlet  helium  are  of  interest.  The 
atomic  state  1 1,0, 0>  Is  negligibly  affected  by  the  elec-  | 


FIG.  1.  Schematic  diagram  for  translHoo.  Wavy  line:  transi- 
tion produced  by  the  radiation  field.  Straight  solid  line:  transi- 
tion produced  by  the  stochastic  electric  field. 

trie  fields,  but  the  n = 3,  J=l,2  states  are  considerab- 
ly affected  by  these  fields. 

For  convenience,  the  states  in  the  >t  = l and  >i  = 3 
level  11,0,0),  13,1,1),  13,1,0),  13, 1,-1),  13,2,2), 
•2>2,1),  13,2,0),  13, 2,- 1),  and  13, 2, - 2)  are  denoted 
by  0, 1, . . , , 8,  respectively. 

If  the  atom  and  radiation  field  are  initially  in  the 
10,  Af|)  states,  the  probability  amplitude  C|(f) 

= C i)(f),  < = 1, , . . , 8 for  finding  the  helium  atom  In 
the  f substate  of  level  n = 3 after  absorbing  a photon 
of  energy  Hui  at  time  / Is  (C,  = C„,,^,) 

^ (S  <«■  \H'\j)  exp«(v,, /)C,(/) 

+ E'  <«l«'l".^>e*P('<*>i.i*»)0C,^„(f)Y 

where  la,Af)  is  now  the  abbreviation  for  \n,J,M,S), 
and  Sts,  »>  ^ fhe  summation  for  all  I a,  N)  except 
< = 0 8. 

Neglecting  far-off  resonant  terms  which  are  small  by 
orders  of  magnitude  compared  to  the  nearly  or  exactly 
resonant  terms  as  well  as  the  terms  associated  with  the 
high-order  optical  transitions  (s  > 8),  and  using  the 
selection  rules  described  by  Eqs.  (24)  and  (25),  one 
can  reduce  Eq.  (26)  to'**** 


= P\  •xrit(u>c  + wt  - a>)flC,(f)  +Ct{[£lsxptf^,)  + <E*,exp(f^J)exp{<(w,  - wi  + 0)1] 

+ [£*, exp(-  <♦.)  + exp(-  ex|>(<(b>,  - u» ^ - n)t])C,(l)  + CfE*{exp(i<i,.) exp(f  (u.,  + 0)1] 

+ exp(-  <♦,)  exp(f (b>,  - 0)f])C,(/)  + Cf{t£j  exp(f^,)  - fF*  exp«0^ J expfKw,  + + 0)f] 

+ [£jsxp(-<^,)-ie*exp(-f*J]exp[f(ai,  +«!-  0)l]}C,(f),  (27) 

= Ft  exp(<(w,-  w)flC,(0  +C|{(£«exp(f(»),)  +fE*  exp(i^J]  exp(<(u>,  - uit  + 0)1] 

+ [E*  sxp(-  <♦,)  + f£} sxp(-  f#^]  »xp[i(u,  - u)t  - 0)f])C,(l)  + Cj£l{exp«^,)  exp[<(w,  + 0)1] 


i 

i 

i 
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+ e*p(-  10.)  exp[i(ui,  - n)/]}C|(0  -•'Cj{[£!exp(<^,)-i£^exp(i4i,)]exp(f'(u|  +0)^  + O)/] 
+ ffijexpl-  «>,)  - i£jexp(- 10,)]  exp(i(u),  + 01^  - n)/I}C,(0, 


= Pi  exp(<(u)o  - <*<1  - u))/lCo(0  + C}{(eJ exp(i4>,)  + i£jexp(i<),)]  exp[i (u),  - Wt  + n)/| 

+ [£!«xp(-  «>,)  + »'£5«*P(“  «*?(<(“!  - «t-  n)'lKi(0  + C5£;{exp(/>,)  explf(u),  + 0)/] 

+ exp(-  i(p,)  exp[/{ta),-  n)/l}C,(/)  + C}{(£j  exp(t>,)-i£jexp(i0,)lexp{/(u>,  + 0)4  + n)/| 

+ (£;  exp(-  - i£}exp(-  i«,)]  exp[r(a(,  + - n)/]}C,(/), 

= - C}*{f£!  exp(r0,)  - i£:Jexp((<>,)l  exp[i (-  ui,  + uii  + «)<]  + [JEj exp(- 10,)  - /Ej exp(-  (>,)] 

X exp(i {-  w,  + oil  - n)/l}C, (/), 

= - C}*£S{exp((<^.)exp[f(-  w,  + {J)r)  + exp(-i<(.,)exp(;(-u),-  n)/l}C,(0 

- Cj*{(£; exp(i>,) - i£;exp(i>^]  exp[<(-  ui,  + ui  + n)(]  + [EjexpC- 1>.)  - i£:Jexp(-  i<t>,)] 
xexpfi(-  ui,  + uij,  - 0)/l}C2(/), 

= - ^1*11^! exp(/<(),)  + <£5exp(i>,)]  oxp[i (-  u>,  - uij,  + n)<|  + [£5  exp(-  i<t>,)  + 1£; exp(-  i((>,)l 
xexp[<(-  u)|  - Ui-  n)/]}Ci(/)-  Cj*£j{exp{i0,)  exp[i(-  a),  + £1)/]  + exp(-  exp(f(-  Wj  - n)£|}C2(0 

- C^*{[£jexp(«^,)-  <£5exp(i>J]exp[j(-  w,  + 02,  + {})/]  + [£5exp(-  «0,)  - i£5exp(- 
X exp{i' (- tai,  + u)i  - n)/]}Cj(/), 

i<'i  - Wi  + fl)/]  + [£j exp(-  10.)  + i£5exp(-  i<t>,)] 

xexp[«(-  o),  - Wi-  n)/]}C2(/)-  C5*£j{exp(i0,)  exp[£(-  w,  + n)<]  + exp(-  t0,)exp{i(-  u)i  - n)<)}Cj(0, 

= - Cj*{[£;exp(/i>>,)  + i£5exp(i>,)]exp(<(-  u),-  wi  + £1)/] 

+ [£;exp(-  + ££;exp(-  i>,)]exp[i(-  oj,  - aii  - £J)£]}Cj(£), 


/>,=  - + (,)i)(e,-ie,)<l|x  + <)i|0)(N/2toT«<.))''’, 

Pt  = - euioe,(2 1 x [ 0)  (Ar/2t,  t/Iw)*  '*, 

P)  = - 4«’(wo  - «>t)(e,  + »e,K3  |x  - iy  \ 0){N/2ttTKu)' 

C}=  (ie/4m  I*  - I'y  1 4),  C}=  (£e/2ff)<l  |x  1 5),  0?=  {ie/mi  |x  + «>  1 6), 

Cj=««/4«)<2|x-iyl5).  C5=«e/2«)<2|x|8),  C}=  (.•e/4«)<2|x +.-y  I?), 

C?  = «e/4«)(3  |x  - I'v  1 6) , Cj  = (£e/2«)<3 1 x 1 7) , Cj  = (fe/4«)<3  |x  + ly  1 8>. 


In  these  equations,  c,  represents  the^  component  of  the  polarization  vector,  C{*  denotes  the  complex  conjugate  of 
C{,  = and  <•>,  = **'(£^,^2)“  ^?a,j))-  Figure  1 Is  a schematic  diagram  showing  the  transitions  corresponding  to 

Eqs.  (27)- (34). 


The  O-'i  atomic  transition  probabilities  per  unit  time  are  given  by 

F,=  |C,(/)|*f-'  tf  = l 8),  (41) 

where  t Is  much  smaller  than  the  lifetime  of  level  0,  but  much  greater  than  2»/£l.  Accordingly,  one  can  write 

Co(0«l,  (42) 

C,(f)«l  (£  = 1 8).  (43) 


It  should  bo  noted  that  even  11  C|(f)«l,  £ = 1 8,  there  can  still  be  an  effect  ct  observable  magnitude  11  a large 

number  of  independent  systems  contribute. 


Substituting  Eq.  (42)  Into  Eq.  (27)  and  integrating  by  parts,  one  obtains 


♦ n.)  * «•  .W(-  CM  * 
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■ti) 


* ■"  '''t((^«p(''0.>  - i-g;  exp(.\»,)  I 

+ [EtexpC-  !♦.)-  /£‘exp(-  " /'  [c|((^. •*?('>.)  + '£*,e*P«^,)l 

+ (£;exp(-.«.)-i£le)q)(- /<.,)!  ^ ] <«•  <♦<) 

Because  o(  the  assumption  on  the  transition  probability  as  described  in  Eq.  (43),  the  second,  third,  and  fourth 
terms  at  the  right-hand  side  of  the  above  equation  are  neglected  in  comparison  to  the  rest.  This  is  also  justifiable 
by  the  direct  estimation  of  these  terms  as  will  be  seen  in  Sec.  IV. 

In  view  of  Eqs.  (30)— (32),  dC^(l)/dt,  dC^U)/dl,  and  dCf{l)/dl  cannot  be  neglected. 

Substituting  Eqs.  (30)— (32)  into  Eq.  (44)  and  differentiating  with  respect  to  I yields 


^-/,(0C,(0=g,«), 

which  has  the  solution 

C,(0  = exp(  /,'/(/')  dl'\ /.‘/r,(/')  exp(  /;>(/')  rf/”  Id/', 
where 


MS) 


M6) 


/i(0=  |Ct|*^[£;exp(i0.)  + if:*,exp(i0,)|[£jexp(i>,)  - /£; eJtp(i((»,)l ^ + i£jexp(-  i<fi.)  + ;£;exp(- 10,)(* 
X [i(u)i  - ta>t  - O)]-'  + I £jexp(«(<>,)  +f£iexp(f(^,)  - Wi  + SJ)]-*  + [£“  exp(-  id),)  + i£jexp(-  (<>,)1 

x(£;exp(-/,J-i£;exp(-/^,)l^^i:;^)  . |cf [»£«; (exp(f2^.)f^  . 

+ exp(-  I exp(id>,)  - /£«  exp(id>^][£i  exp(j>,)  + j£;  exp(id>,)] 

+ |£jexp(-i(t>,)-  i£5exp(-/«,)|’[»(w,  + (1)^-  0)]’*+  |£5exp(i(<),)  - t£' exp(f</),)|’[i((Oi  + + 0)]*' 

+ f£;exp(-  f0,)-  ,£‘,exp(-  ,>,)K£jexp(-  id,.)  + f£;exp(- 


M7) 


and 


<f,(/)  = exp(i(uj  + u)i  - w)/]  + ^C}C}*£»  ^exp(»((>.)  + exp(- 

x{[£jexp(t<»,)  - i£*  exp(i>,))  exp(i(-  u,  + Wi  + fl)/]  + [£jexp(-  id.)  - i£5exp(-  i«,)]  exp[i(-  i),  + Wi  - fl)/]} 
.C?C!-£;((£;exp(f,.)-  i£iexp(/,Jl^ieji^^  4 [£;exp(-  .>.)-  .•£;exp(- .-,,), 

x{e:q)(i>,)  exp(«(-  u,  + fl)/]  + exp(-  id>,)  exp(«(-  u»,  - n)/]}jc,(f)  + C}Cf*^(£5  exp{id>.)  - i£}  exp(i<f ,)] 
x{(£jexp(«^,)  - f£5exp(i>,)]  exp(j(-  u>,  + + tl)/]  + (£}e)q>(-  id>,)  - f£jexp(-  i«,)|  exp(«(-  u),  + ui^  - n)/)jC,(f). 


(48) 


Because  of  the  random  properties  of  phase,  the  direction,  and  the  amplitude,  and  the  high  frequency  (Hf » 1)  of 
the  stochastic  electric  fields,  J^/iU')dt'  fluctuates  negligibly  about 


(49) 


where 


' 3^  V 2[(u(|  - uii)’  - fl*)  u)}  - n*  2[((»), + (iii)’- n^1  / 


38S5 
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(50) 


r ' 

i 


)■  ■ 


= / — 6(wi-  toiz.)  , . <^i  * _\ 

ii*  \ (ufi  - — 0*  w}  - 0*  + u>|,)*  - 0*/  * 


In  these  equations*  ( > means  the  average  over  the  phase*  the  direction,  and  the  amplitude  oC  the  stochasUc  fields* 
denotes  the  Bohr  radius*  Q = 6. 7,  and  ofi  Is  204  cm**  in  ware-number  units. 


Accordingly*  one  obtains 

C, (0-exp(-<Aui,0/,'#i(nexptfAu), /')<«' •/>i«xp(- 1 Aui,0  //exp(<(ai,  + <.)j,  + Au),  - 


and 


T]  = 2»|i’i  + Uli  + Aui|  - w). 


In  theae  equations,  the  second  and  third  terms  o(  Eq.  (48)  are  neglected  compared  to  the  first  term.  It  will  be 
Justified  In  Sec.  IV. 


(51) 

(52) 


Equation  (52)  means  that  the  0 — 1 allowed  line  (a,  line)  is  shifted  as  much  as  Aw,,  which  is  proportional  to  the 
mean-square  stochastic  electric  field  (the  quadratic  Stark  effect).  Similarly,  the  probability  amplitudes  and  the 
transition  probabilities  for  the  2 and  3 levels  (w  and  o.  lines)  are 


(53) 

(54) 


(55) 


(56) 

(57) 


(58) 


C|(0  = P}  exp(-  fA<i)]()  exp[i(u)|  + Au>] - ui)/']df', 
r|  = 2*  |P]  |*9((i)|  + AW]- w), 
where 

* 2[(u,, -«,)*- n’l  ^ u.!-n»  ^ 2[(u.;+wi)‘-n*l  j 

_ ^e*(E*)ai  / 3((4>i  - (i>t.)  4a)i  3(u)i  wt)  \ 

36*  \ (u>i  -u>i)’  - n*  wj-n’  (u>,  + (lij)*  - n* j ’ 

Ci(<)  = PjSap(-  iA(i>]<)  exp(f  (wg  - ui  + au)|- 

rj  = 2»|/’||*9(wj-  + Alalg-  b)), 

where 

* W \ 2t(b.,-u,i)‘-n»]  b.{-n‘  * 2(u.,  + wj«-n»  J 

= + 3aL  4 N 

3i*  \ (a»|  - - n*  u>}-  0*  (wj + 0?^)*- O*/ ' 

Substitution  04  Eq.  (51)  into  Eq.  (30)  yields 

[(*!•■*««  - 

X /e«p[<(a>.-b)i-t-2h)i.-n-h))<]-l  , e)q?[»(-h»,-t-u)i-Aui.-ni<1-l\-| 

V Uig  - bi| -f  2(0^- n- w - uii Awi  - n /J  ' 

Accordingly,  the  probability  for  the  transition  0 — 4 Is 

" ^W((i)*  +il*  + - u),  + 2b)i  + n - w)  + 8(u»g  - 01,  + 2wi  - n - ui)]  (-O»i  + O)i-  Ao),±n#0). 

similarly  one  obtains  for  the  probabilities  for  the  transitions  0 — 5,  0 — 6,  0 — 7,  0 — Bas 

- 6)P  A(o.g  + o.t  + Ao),  -o.)*  ~2(u)'g  V A0),-  0>)»  '<"» - u),  + 0>i  + fJ - o,)  + 9(a.g - o.,  + u.^  - 0 - o.)],  (61) 

2^|^/|Ptl«|<ll«-t-fyl6>|»  ^ IP.|M<8I»-<vI6>I»  IP«l»l(2lxl6)l»\ 

6S*  V 2(u)g  + 01^  + Ahl|  - U>)’  2(h)g  - 0»1  + AUg  - (Jp  (h>g  + Ah)]  - hi)*  / 

*(6(hlg- h)j  + 0- h»)  + 6(hlg- h)i  - 0- hi)),  (**) 

**  l^dg'?.!,*  [9(h»g-h»,-2h»i  + n-h))  + 6(h)g-h),-2h.i-n-h))l 


(59) 


(60) 


(64) 
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In  these  equations,  we  used  |P||’^  IPjl’ = which 
comes  from  the  fact  that  the  radiation  Intensity  does  not 
change  significantly  In  the  range  of  u>i  2ti>i  about  uig. 

In  deriving  Eq.  (62),  we  neglect  P^Pf  and  PfPj  since 
the  average  of  the  cross  terms  between  the  different 
frequency  components  with  no  phase  relations  of  the 
radiation  field  is  zero. 

Equations  (60)— (64)  mean  that  there  are  two  satellites 
disposed  symmetrically  In  pairs  about  each  of  the 
forbidden  lines,  0-4,  0-5,  0-6,  0-7,  and  0-8, 
and  separated  from  them  by  (1. 

It  Is  Important  to  point  out  that  the  spectral  shift  of 
the  satellite  line  Is  not  seen  In  Eqs.  (60)— (64)  ("asym- 
metrical Stark  repulsion”).  The  spectral  shift  of  the 
satellite  line  can  be  seen  in  the  case  of  the  selective 
excitation  of  the  forbidden  transition  by  the  laser.  This 
result  cannot  be  explained  by  the  coupled  equations  with 
the  phenomenological  damping  constants  of  Prosnltz 
et  al. ' Similar  results  have  been  found  by  Grlem‘*  since 
Grlem  has  not  treated  the  laser  radiation. 

The  ratios  of  the  absorption  coefficients  at  these 
satellite  frequencies  to  that  of  the  allowed  line  are 


+ AW]  ± fl)* 


(65) 


51^  »e^aS<E*)  / 1 ^ 1 \ 

* W*  \ (uij  + ± Cl)’  (u>i  ~ + 1 n)T  / ' 
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fo.  <{,  0 — 5,  0 — 6,  0 — 7,  and  0 — 8 satellites,  re- 
spectl  'fly,  vhere  the  minus  sign  In  the  above  formulas 
corresp  to  the  near  satellite  and  the  plus  sign  to  the 
far  satellites. 

ff  B = 0 (u)i  = 0),  there  are  only  two  satellites  dis- 
posed symmetrically  in  pairs  about  a forbidden 
1 'S,-3  'Oj  line  and  separated  from  It  by  0.  The  ratio 
Sg  defined  above  In  this  case  Is 

„ 5gg*a|(g*> 

«*(ai7+^±  n)»  ’ 

where 


(«i.t  = W = 67) 


02) 


has  a factor  and  a denominator  which  are  different  from 
Eq.  (71)  by  1,  and  by  nwvO,  respectively. 

The  perturbation  calculations  by  Baranger  and  Mozer' 
and  Grlem*’  could  not  include  the  detailed  coupling 
mechanism  between  the  upper  states  of  the  allowed  and 
forbidden  transitions  and  the  Stark  shifts  of  the  allowed 
lines.  The  deficiency  In  the  coupling  mechanism  gives 
rise  to  the  above  difference  in  factor  and  the  fault  that 
does  not  include  the  Stark  shift  to  Au  v 0 of  the  denomi- 
nator. The  three- level  rotating- wave  approximation  of 
Prosnltz  el  al.‘  has  the  former  deficiency  as  the  pertur- 
bation theory  and  would  give  rise  to  a different  denomi- 
nator since  their  theory  would  claim  Stark  shifts  of  the 
satellites  in  the  absence  of  the  laser  radiation. 


III.  CONCLUSIONS 

A theory  of  the  spectral  shifts  of  o„  a.,  and  v com- 
ponents of  the  allowed  line  (1  ‘So-3'P,)  and  the  satellites 
about  the  forbidden  line  (1  ‘Sg-S  'Hj)  of  a He  I plasma  In 
the  static  magnetic  field  has  been  developed  from  first 
principles.  In  which  the  atomic  states'are  exposed  to 
(longitudinal)  turbulent  electric  fields  and  (transverse) 
radiation  fields.  The  o.,  o.,  and  n components  of  the 
allowed  line  are  shifted  as  much  as  Au,,  Auj,  and  Auj, 
respectively,  while  the  sateUlte  lines  are  not.  The 
ratio  of  the  Integrated  intensity  of  the  satellite  line  to 
that  of  the  allowed  line  is  different  from  the  ratio  ob- 
tained from  the  second-order  perturbation  theory*  or 
the  three- level  rotating-wave  approximation.  * From  the 
observation  of  the  distance  in  frequency  between  two 
satellites  of  a forbidden  line  and  of  the  shifts  of  the 
allowed  lines  or  the  ratio  of  the  integrated  Intensity  of 
the  satellite  line  to  that  of  allowed  line,  a quantitative 
determination  of  the  frequency  and  intensity  of  the 
stochastic  electric  field  Is  possible. 

IV.  VALIDITY  OF  THEORY 

We  discuss  as  to  whether  the  solutions  for  C,(0, 

Cj(/),  and  Cj(f)  given  In  Eqs.  (51),  (53),  and  (56), 
respectively,  are  good  approximate  solutions  of  the 
coupled  equations  (27)- (24).  In  the  following,  we  denote 
Cg(0,  I = 1,2, 3 given  In  Eqs.  (51),  (53),  and  (56)  as 
Cf(0,  f=  1,2,3  In  order  to  distinguish  from  the  exact 
solutions  C,((),  I = 1,2, 3. 

(1)  First,  we  prove  that  CJ(I)  satlslles  the  following 
approximate  equation: 

^§^-/t(0C?(f)-g,*(0,  (73) 

where  g|*(0  can  be  expressed  by  replacing  C|(0  with 
CJd)  In  Eq.  (48). 


It  should  be  pointed  out  that  S.  given  by  Baranger  and 
Mozer*  as 


Since  C^(()  satisfies 

= P^  exp(«(ta)g  + fi,  - w)/],  (74) 

the  validity  of  Eq.  (73)  depends  on  that 
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x{(£jexp(/<>,)-  <ff}e*ptf^Jlexp(<(-  <•>,  + <i»i  + n)<'|  + (£*exi)(-  »>,)-  /£'*exp(-*-  /(p,)lexp(i(-  u»,  + wt-  n)/'|} 
-*•  <1Cl*E«^[E« exp(>0,) - lE* exp(<»  J| ^ * l^®****"  '♦.)-  '^,e*P<-  «>,)l 
* ^•*P*'*^«*  •»?(<(-  <*>i  + n)/'l  + «xp(-  !♦,)  exp{((-  u),  - n)/'lljcj(/') 

J + [gJexpC-  iip,)-  <£^exp(-i<>,)j 

^ + ■ igtexp(i>Jlexp(i(-  u,  + uij,  + tl)/') 

+ [£jexp(-  !>,)-  i£jexp(-  i^,)]exp(«(-  w,  + u)j,-  n)/'l}c;(/')|exp(i  A<i),/')<</' j 
can  be  Ignored  In  eompartaon  to 

I Cf  (0 1 * = I /o‘  J®,  expfi  (w,  + Ui  + A(j,  - (i»)/)  d/ 1 * = 2W I Pi  I *, 


(75) 


(76) 


where  the  summation  with  reapect  to  w has  been  made  in  Eq.  (76),  since  the  radiation  field  covers  a wide  spread  of 
(requenclss  and  the  contributions  from  various  frequencies  are  additive. 

From  the  direct  calculation  of  Eq.  (75)  alter  substituting  C{(f)  and  C^(f), 

\C\\\  1C!1‘,  |C1|»,  lC}|‘MQe*ai/>f‘) 
and 

(ill,  (1,  ^1, 

for  most  experimental  cases,  we  can  be  sure  that  if 


1 1 1 

((ill,  + 20]  + 0(i)|  - A(ii})'  ’ ((Vi  - 20  + A<iii  - A(i>])’  ’ (dll  + ^i  - Ouij)*  (Qo5)V‘<eV 
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and 


1 1 1 

(2(1)  I + Auii-  Odij)*  ’ (2uii  + 20  + &(ii|  - A(i)|)*  ’ (2(iii  - 20  + O&ii  - A(i)j)’  ^ (Qao)’*A(E’)*>  ’ 


(80) 


ki(f)  can  be  ignored  in  comparison  with  ICf(f)l*.  Inequalities  (79)  and  (80)  are  the  same  as  the  criteria  for  the 
validity  of  the  perturbation  calculation,  which  can  be  stated  in  a crude  form  that  and  I8uii  - 280 1 are  greater 
than  the  perturbing  energy  of  the  electric  field  £,,  where 

g»~(109)''’«o((E*»‘'* 


(li)  We  proved  that  C\{t)  satisfies  Eq.  (73).  Hence,  C)(f)  obviously  satisfies 

where 


(81) 
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V{t)  = - J ^C}  ^[E^exp«<»,)  -t-  «E*  exp(t» J]  + [g*exp(- tifi,) +fE*exp(- 10,)] 

+ Ct^[E; exp(f »,)  - fg*  exp«» J]  " «^te*P(- 

„exp{f((v.-t-wi-0)<lWC|(f)1j 

f((.i,+(Vi-0)  ) dt  J ’ 


(83) 


where  dC\{t)/dt,  dC\{t)/dt,  and  dCt(t}/dt  is  given  by  replacing  C,(/),  i = l,2,3  by  CJ(/),  i = l,2,3  in  Eqs.  {30)-(S2), 
respectively. 
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Accordingly,  if 

|q(0l’-lv(/)l». 


V{,) f,((£;exp(t>.)  * ^ [£.  gaL«^^irn_-«v|)  < 

X CfEj  ^expdd.,)  + e;q,{-  + q((£;exp(.>.)  - ,£;exp(,0,)) 


, exp(;(u)i  +u)L  + a)l 


that  Is,  C^U)  Is  an  approximate  solution  of  Eq.  (44). 

From 

lq(/)|’,  |Cf(ol’,  (8' 

and  Eq.  (77),  we  can  be  sure  that  If 

then  Eq.  (84)  Is  valid,  that  Is,  C}(f)  Is  a good  approximate  solution  of  Eq,  (27).  Similarly,  we  can  prove  the  same 
for  the  remaining  equations. 
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